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ABSTRACT 


This  thesis  lays  the  foundation  for  the  dynamic  stall 
Investigation  being  conducted  at  the  Fluid  Mechanics  Lab- 
oratory at  NASA-Ames  Research  Center.  Using  existing  optical 
and  electrical  equipment,  a  new  dedicated  Micro-VAX  computer 
and  Labstar  software,  an  indraft  transonic  wind  tunnel  and 
able  technicians  to  make  the  proper  Interface  hardware,  the 
project  came  together  in  a  new  test  facility  at  the  Fluid 
Mechanics  Laboratory.  The  goal  of  the  thesis  was  to  obtain 
both  qualitative  and  quantitative  information  about  the  wake 
profiles  of  an  airfoil  in  steady  state  operations  at  varying 
angles  of  attack  and  tunnel  conditions.  To  accomplish  this 
task,  schlieren  photography  was  used  to  obtain  a  qualitative 
picture  of  the  flow  field.  With  this  Information,  a  two  com- 
ponent Laser  Doppler  Velocimeter  was  set  up  to  accurately 
measure  the  velocity  profiles  that  correspond  to  the  schlie- 
ren photographs.  Once  this  preliminary  work  is  completed, 
the  same  apparatus  will  be  used  to  further  Investigate  the 
unsteady  dynamic  stall  phenomenon. 
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I.  INTRODUCTION 

A.   DESCRIPTION  OF  THE  OVERALL  DYNAMIC  STALL  PROJECT 

Static  stall  Is  a  well  understood  and  documented  phenom- 
enon. A  wing  or  rotor  blade  that  is  operating  at  a  high  a, 
angle  of  attack,  will  experience  a  point  where  separated 
flow  will  cause  an  accelerated  decrease  in  lift  and  a  large 
rapid  increase  in  drag.  This  a  at  which  the  wing  will  pro- 
duce the  maximum  lift,  Cjmax,  ^3  a  very  important  aspect 
of  airfoil  performance.  This  stall  a  determines  the  stall 
speed  of  the  aircraft.  Considerable  effort  has  been  expended 
to  find  ways  to  decrease  the  stall  speed  by  Increasing  the 
Cjmax.  The  majority  of  this  effort  has  been  directed  to- 
ward a  gradual  approach  to  the  stall  a,  in  other  words, 
static  stall  research. 

The  wing  can  be  made  to  continue  to  produce  lift  past 
the  static  stall  a  if  the  wing  angle  of  attack  Is  in- 
creased rapidly  past  the  stall  a.  The  phenomenon  of 
dynamic  stall  is  manifested  in  the  rapid  decrease  of  the 
lift  at  the  higher  stall  a.  There  are  four  basic  stages  of 
dynamic  stall  that  have  been  documented  [Ref.  13: 

1.  Tral 1 lng-edge   flow  reversal,  progressing  forward  along 
the  alrf ol 1 . 

2.  Formation,  growth  and  shedding  of  the  vortex. 

3.  Ful 1  stal 1  . 

4.  Boundary  layer  reattachment. 

The  understanding  of  this  phenomenon   will  be  beneficial  to 
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the  helicopter,  where  constantly  oscillating  rotor  blades 
are  consistently  in  the  dynamic  stall  region  throughout 
their  revolution,  and  jet  fighters,  where  agile  and  rapid 
maneuvering  through  the  stall  a  is  required  for  surviva- 
bility. By  quantifying  the  velocity  field  around  and  behind 
the  airfoil,  a  better  understanding  of  the  dynamics  of  this 
phenomenon  will  be  gained.  This  new  analytical  information 
will  provide  valuable  Information  to  the  design  engineer  to 
expand  the  flight  envelope  of  future  fighter  aircraft  and 
improve  the  performance  of  the  helicopter  rotor  system. 

B.   THESIS  GOALS 

This  thesis  lays  the  foundation  for  the  dynamic  stall 
investigation  project  being  conducted  at  the  Fluid  Mechanics 
Laboratory  at  NASA-Ames  Research  Center.  Using  existing 
optical  and  electrical  equipment,  a  new  dedicated  Micro-VAX 
computer,  new  Labstar  software,  an  existing  transonic  in- 
draft wind  tunnel  and  able  technicians  to  make  the  proper 
interfaces,  the  project  came  together  in  a  new  test  cell  of 
the  Fluid  Mechanics  Laboratory.  The  goal  of  this  thesis  was 
to  first  obtain  stroboscopic  schlieren  photographs  of  a  sta- 
tionary airfoil  at  different  a's  to  obtain  qualitative 
data  on  the  flow  past  the  airfoil  and  the  size,  shape  and 
extent  of  the  wake  and  then  set  up  a  Laser  Doppler  Veloci- 
meter  to  obtain  quantitative  velocity  data  in  order  to  con- 
struct an  accurate  velocity  profile  of  the  flow  field. 
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C.   SCHLIEREN  AND  LASER  DOPPLER  VELOCIMETRY  BASICS 

1 .   Schl leren 

Many  optical  flow  visualization  techniques  have  been 
developed  for  gas  flows.  Each  of  these  techniques  are  sensi- 
tive to  variations  in  density  and  derivatives  of  density  In 
the  flow.  Interf erometry  Is  sensitive  to  the  density  of  the 
gas  flow;  f.  Schl leren  is  sensitive  to  the  density  grad- 
ient or  the  change  in  density  with  respect  to  distance; 
SfSy.  Shadowgraphy  is  sensitive  to  the  second  derivative 
of  the  density;  6^^^iy^.  All  the  techniques  rely  on 
the  fact  that  the  speed  of  light  is  a  function  of  the  local 
index  of  refraction  in  that  medium.  In  gases,  the  index  of 
refraction  is  a  function  of  density. 

The  schl leren,  German  for  streak,  technique  allows 
one  to  view,  qualitatively,  the  density  gradient  in  a  given 
flow  field.  This  visualization  system  is  considerably  sim- 
pler and  cheaper  to  use  than  other  techniques.  Therefore,  it 
was  chosen  as  the  preliminary  technique  to  obtain  qualita- 
tive information  about  the  air  flow  around  and  in  the  wake 
of  the  airfoil.  The  schlieren  setup  involves  using  a  line 
source  of  light  and  projecting  it  through  the  test  section 
where  the  flow  field  alters  the  path  of  the  light.  The  line 
source  Is  usually  a  monochromatic  gas  discharge  tube.  The 
light  is  projected  through  a  narrow  slit  to  provide  a  line 
of  light  of  uniform  illumination.  The  dimensions  of  the  slit 
control  the  definition  of   the   Image   seen  on  the  screen.  A 
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high  quality  concave  mirror  Is  usually  used  to  direct  the 
light  from  the  source  through  the  test  section.  On  the  other 
side  of  the  test  section,  another  similar  mirror  Is  used  to 
reflect  the  light  to  a  viewing  screen  or  camera.  Each  point 
of  the  line  source  will  illuminate  each  point  in  the  test 
section  if  the  source  Is  col  1 Imated  and  of  equal  optical 
distance  on  both  sides  of  the  test  section.  The  light  from 
the  second  mirror  is  focused  onto  a  knife  edge.  The  place- 
ment and  adjustment  of  the  knife  edge  Is  very  critical  as 
the  sensitivity  of  the  entire  system  depends  on  It.  The 
knife  edge  is  placed  In  the  plane  of  the  source  Image  and 
interferes  with  a  portion  of  this  image.  A  portion  of  the 
light  from  each  point  in  the  test  section  is  blocked  which 
decreases  the  light  impinging  on  the  viewing  screen.  The 
density  gradients  encountered  in  the  test  section  normal  to 
the  knife  edge  cause  a  refraction  of  light  allowing  more 
light  from  that  portion  of  the  test  section  to  pass  over  the 
knife  edge  and  Illuminate  the  viewing  screen.  By  rotating 
both  the  line  source  and  the  knife  edge  both  directional 
derivatives  of  density  can  be  viewed.  Also,  through  the 
lateral  adjustment  of  this  knife  edge  at  the  source  image, 
a  fine  control  over  the  degree  of  detail  resulting  from  this 
technique  can  be  exercised. 

The  viewing  screen  can  be  any  image  recording 
device.  Polaroid  cameras  can  be  set  up  for  Instantaneous 
results  or  35  mm  cameras  can   be  mounted  at  the  focal  length 
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of   the   lens   to  obtain   numerous  photographs  during  the 
testing.  Video  cameras  can  also  be  used. 
2.   Laser  Doppler  Veloclmetrv 

The  techniques  for  measuring  velocity  in  the  past 
used  sensors,  which  need  to  be  Introduced  into  the  flow. 
This  necessarily  meant  Intruding  into  and  affecting  the 
flow.  During  the  last  two  decades,  a  new  technique,  Laser 
Doppler  Veloclmetry  or  LDV,  has  been  available.  LDV's  are 
Indeed  a  better  technique.  It  is  nonlntruslve  and  is  capable 
of  directly  measuring  velocities  of  particles  In  the  flow. 
Laser  Doppler  Veloclmetry  has  been  very  successful  In  meas- 
uring flow  velocities  in  turbulent  and  reverse  flows.  It  has 
been  successfully  used  for  the  past  two  decades  to  obtain 
velocity  measurements  accurately. 

Laser  Doppler  Veloclmetry  uses  two  coherent*  light 
beams  from  the  same  source,  which  intersect  to  form  the 
measurement  volume.  The  two  light  beams  form  an  Interference 
fringe  system  that  is  crossed  by  particles  moving  with  the 
flow.  The  light  scattered  by  the  particles  Is  Doppler 
shifted2.  The  frequency  of  this  scattered  light  is 
directly   proportional   to    the    particle   velocity.   With 


1.  Coherent  Light  =  light  in  which  the  electromagnetic 
waves  maintain  a  fixed  phase  relationship.  (Random  House 
College  Dictionary) 

2.  Doppler  Shift  =  the  magnitude  of  the  apparent  change 
In  frequency  resulting  from  the  relative  motion  of  the 
light   source   and  the  receiver.  (Random  House  College  Diet) 
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suitable  selection  of  flow  particles,  known  also  as  seeds, 
the  fluid  velocities  can  be  measured  from  the  particle 
velocities.  Using  photodetectors,  the  light  scattered  by  the 
particles  is  converted  into  electrical  signals.  This  signal, 
after  processing,  provides  frequency  information  from  which 
the  velocity  can  be  determined.  The  measurement  of  two 
components  of  the  velocity  can  be  obtained  by  setting  up  two 
orthogonal  systems  with  laser  beams  of  two  colors,  In  the 
same  optics  train,  thus  Identifying  each  component  with  a 
particular  color. 
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II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

A.   SCHLIEREN  SETUP 

The  achlieren  system  for  this  experiment  was  set  up 
using  both  a  continuous  light  source  and  a  stroboscoplc 
source.  The  continuous  source  was  used  to  provide  a  time 
averaged  schlleren  photograph  and  also  to  allow  use  of  a 
video  tape  recorder  to  record  the  time  history  of  the  flow. 
To  provide  Instantaneous  schlleren  photographs  of  the  flow 
field,  the  stroboscope  was  used.  See  Figure  2.1  for  the 
schematic  of  the  setup. 

To  obtain  the  35  mm  photographs  of  the  flow  field,  a  US 
Scientific  Strobrlte  Stroboscope  was  used.  The  strobe  had 
five  settings,  20  Hz  to  1  KHz,  to  adjust  the  duration  of  the 
strobe.  The  70  Hz  setting  proved  to  be  the  optimal  setting. 
A  mechanically  adjustable  silt  In  front  of  the  strobe  was 
used  to  produce  a  line  source  of  light  needed  for  the  tech- 
nique to  operate  effectively.  A  large  concave  mirror  with  a 
305  cm  focal  length  reflects  this  light  sheet  towards  the 
test  section.  The  test  section  had  optical  access  through 
both  sides  of  the  tunnel.  The  test  section's  dimensions  were 
13  cm  by  25  cm  and  was  constructed  out  of  2.5  cm  thick 
plexlglas.  The  light  passing  through  the  test  section  was 
then  reflected  by  an  identical  mirror  located  opposite  the 
first   mirror  on  the  other  side  of  the  test  section  as  shown 
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Top     View 


Concave  Mirror 
(120"  Focal  Length) 


Vertical  Slit 


Stroboscopic 
Light  Source 


Airflow 


Concave  Mirror 
(120"  Focal  Length) 


+35mm  Camera 
■  200mm  Lens 

^  ^Vertical  Knife 
Edge 


Figure  2.1   Schlleren  Schematic 

In  Figure  2.1.  Light  from  the  second  mirror  was  focussed  on 
a  vertical  knife  edge  whose  position  could  be  adjusted  to 
get  a  highly  sensitive  system.  The  knife  edge  was  oriented 
vertically  to  obtain  the  visual  signature  of  the  density 
gradient  In  the  vertical  plane  of  the   test   section.   After 
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the  knife  edge,  a  200   mm   lens  was. used  to  focus  the  image 
onto  the  photographic  device. 

Since  the  test  section  is  25  cm  wide,  the  photographs 
obtained  from  the  system  actually  represent  density  gradient 
effects  averaged  over  the  test  section  width,  as  is  always 
the  case  with  the  schlieren  technique. 

A  Nikon  35  mm  camera  body  with  motordrive  (but  no  lens) 
was  used  to  record  the  schlieren  images.  The  camera  was 
placed  at  the  focal  point  of  the  200  mm  lens.  Fine  adjust- 
ments to  focus  the  image  in  the  camera  were  made  by  a  micro- 
meter slide  on  the  camera  mounting  mechanism.  Pictures  were 
taken  by  opening  the  shutter  with  a  remote  shutter  release 
cable,  before  the  strobe  was  fired  remotely  and  then  closing 
the  shutter.  With  a  little  practice,  this  procedure  was 
quite  rapid.  Photographs  were  taken  at  two  different  Mach 
numbers,  0.30  and  0.45,  and  at  three  different  a's;  0°, 
7.7°,  15.5°.  See  Figures  A.l  through  A. 6,  Appendix  A, 
for  examples  of  schlieren  photographs.  This  series  of  photo- 
graphs was  studied  to  analyze  the  flow  qualitatively  and  to 
determine  the  locations  in  the  wake  for  obtaining  the  veloc- 
ity data  for  the  two  different  Mach  numbers. 

B.   LASER  DOPPLER  VELOCIMETER  SETUP 

This  section  will  discuss  the  specific  LDV  system  setup 
In  test  cell  number  four  at  the  Fluid  Mechanics  Laboratory. 
The   following   sections  will   be   discussed:   transmitting 
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optics,  receiving  optics  and  photodetectors,  test  section, 
signal  processing  and  data  processing.  Figure  2.2  shows  a 
photograph  of  the  transmitting  optics  and  Figure  2.3  shows  a 
photograph  of  the  receiving  side.  Individual  schematics  will 
be  discussed  in  detail  within  the  following  sections. 

The  laser  used  Is  a  Spectra  Physics,  Model  164,  5  Watt 
Argon-Ion  Laser  capable  of  outputting  light  from  0.351  to 
0.528  microns  and  up  to  twelve  discrete  wavelengths  simul- 
taneously CRef.  2:p.  233.  By  choosing  the  blue  beam,  0.4880 
microns,  and  the  green  beam,  0.5145  microns,  for  the  two 
component  Laser  Doppler  Veloclmetry  System,  nearly  75%  of 
the  Initial  output  power  could  be  used.  This  standard 
practice  was  adapted  in  this  case. 

1.   Transmitting  Optics 

Owing  to  traversing  requirements,  the  laser  was 
higher  on  the  optics  table  than  the  10.8  cm  optical  height 
of  the  transmitting  optics.  A  set  of  two  mirrors  at  45° 
angles  to  the  direction  of  the  incident  light  were  used  to 
lower  the  beam  to  the  require  height.  Figure  2.4  shows  a  de- 
tailed schematic  of  the  entire  transmitting  optics  system. 

The  vertically  polarized  beam  exiting  the  multiline 
laser,  after  being  directed  by  mirrors,  Is  passed  through  a 
beam  collimator,  part  of  the  TSI  Model  9106  set.  Incident 
light  is  col  1 Imated  if  the  beam  diameter  does  not  change  as 
It  passes  through  space.  In  reality,  however,  it  Is  suffic- 
ient  if   the  collimator  can  insure  this  quality  through  the 
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Figure  2.2  Transmitting  Optics  Photograph 


Figure  2.3   Receiving  Optics  Photograph 
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Figure  2.4  Transmitting  Optics  Schematic 

relatively  short  optical  path  of  the  system.  Proper  col  11- 
mation  Is  said  to  be  achieved  when  the  waist1  of  the  beam 
coincides  with  the  focal  point  of  the  transmitting  lens.  By 
sliding  the  collimator  sleeve  in  and  out  and  observing  the 
beam  "spots"  on  a  distant  wall,  the  collimator  can  be  set. 


1.  Waist  =  the  point  on  a  beam  where   the   diameter   Is 
the  smal lest . 
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From  the  collimator,  the  beam  passes  through  a  po- 
larization rotator  and  an  attenuator.  The  polarization  ro- 
tator changes  the  polarization  of  the  Incoming  laser  beam  to 
conform  to  the  polarization  requirements  of  the  next  optical 
component,  a  dispersion  prism.  The  attenuator  works  in 
tandem  with  the  rotator  to  insure  passage  of  only  the  horiz- 
ontally polarized  light  through  these  components.  The  dis- 
persion prism  separates  the  wavelengths,  or  colors,  of  the 
light.  The  two  primary  colors  of  this  laser,  green  (0.5145 
microns)  and  blue  (0.4880  microns),  are  used  to  form  the  two 
component  system.  A  set  of  mirrors,  TSI  Model  9107,  is  used 
to  direct  each  of  the  beams  to  the  appropriate  axis  for  the 
transmitting  optics  as  shown  in  Figure  2.4.  As  can  be  seen 
from  the  figure,  the  green  beam  Is  directed  down  the  center- 
line  of  the  optics  and  the  blue  is  directed  25  mm  to  the 
right  of  the  green. 

The  centerline,  green,  beam  passes  through  a  polar- 
ization rotator,  TSI  Model  9102,  to  rotate  the  beam  polari- 
zation 90°,  then  a  beam  splitter,  TSI  Model  9115.  The 
polarization  change  is  necessary  to  Insure  proper  splitting 
into  two  beams  from  the  incident  light.  The  beam  splitter 
divides  the  incident  beam  into  two  beams  of  equal  Intensity 
and  parallel  to  the  optical  axis,  spaced  50  mm  apart. 
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Since  the  flow  field  to  be  Investigated  will  have  a 
turbulent  wake  with  possible  separated  flow  regions,  a  freq- 
uency shift  system  is  used  for  each  of  the  two  colored 
beams.  Upshlftlng  of  the  frequency  of  the  laser  light  prior 
to  entering  the  test  section  will  allow  detection  of  the 
direction  the  particle  traverses  through  the  probe  volume. 
The  shift  results  In  a  moving  fringe  pattern  that  will  gen- 
erate a  photodetector  signal  equal  to  the  shift  frequency 
when  measuring  a  particle  with  zero  velocity.  This  system 
allows  a  greater  range  of  measurements  of  negative  veloc- 
ities and,  thus,  is  required  for  this  system  to  obtain  valid 
measurements. 

The  pair  of  green  beams  now  enter  a  frequency  shift- 
ing device  called  the  Bragg  cell.  This  device  is  an  acousto- 
optlc  unit  that  upshifts  the  frequency  of  an  incident  beam. 
The  shift  is  set  at  40  MHz.  This  device  Is  placed  In  the 
optical  path  of  one  of  the  two  beams  forming  the  probe  vol- 
ume for  that  color.  The  outgoing  beam  will  be  shifted  by  a 
40  MHz  shift  relative  to  the  other  beam.  The  Bragg  cell 
operates  by  deflecting  one  beam  In  a  fan-like  pattern  in 
the  plane  Intersecting  both  Incident  beams.  Each  separate 
deflected  beam  has  a  frequency  that  Is  a  multiple  of  40  MHz 
from  the  Incident  beam.  The  Bragg  cell  should  be  adjusted 
to  shift  the  light  Intensity  to  the  first  order  beam.  This 
Is  accomplished  by  tilting  the  Bragg  cell  In  its  mount. 
This   first  order  beam  should  have  an  intensity  greater  than 
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90%  of  the  Incident  light  beam  [Ref.  3:p.  1483.  An  optical 
wedge  is  Included  in  this  device  at  the  output  end  of  the 
Bragg  Cell  to  steer  the  deflected,  first  order  beam  parallel 
to  the  unshlfted  beam.  The  result  of  this  Bragg  shifting  is 
an  interference  pattern  in  the  probe  volume  that  will  move 
relative  to  a  fixed  particle.  This  moving  fringe  pattern  can 
be  used  to  differentiate  the  direction  of  motion  through  the 
probe  volume.  Particles  moving  with  the  direction  of  motion 
of  the  fringes  gives  a  lower  frequency  and  against  the 
fringe  movement  gives  a  higher  frequency  as  shown  in  the 
following  equation.  The  particle  moving  through  the  probe 
volume  will  produce  a  signal  whose  frequency  Is 

fslgnal  =  '^shift  +/~  fDoppler' 

where  the  sign  is  +  for  particles  moving  in  the  opposite 
direction  to  the  fringe  movement  [Ref.  2:p.  673. 

For  optic  tube  support  purposes  a  TSI,  Model  9176, 
Support  is  used.  This  device  also  serves  to  Increase  the 
effectiveness  of  the  beam  stops  of  the  blue  beam  by  inc- 
reasing the  deflected  beam  separation  distance. 

The  off-axis  blue  beam  Is  now  brought  to  the  center- 
line  using  a  beam  dlsplacer,  TSI  Model  9174.  This  model  em- 
ploys a  rhomboid  prism.  It  displaces  the  beam  a  fixed,  25 
mm,  distance  to  the  centerllne  of  the  optics  while  keeping 
the  beam  parallel  to  the  axis.  The   blue   beam   then   passes 


24 


through  the  same  series  of  optical  components  as  the  green 
beams:  polarization  rotator,  beam  splitter  and  Bragg  cell. 

The  next  two  elements  In  the  optics  train  are  beam 
steering  modules,  TSI  Model  9175,  one  for  each  of  the 
colors.  These  modules  allow  for  steering  of  the  shifted  beam 
to  Insure  proper  beams  crossing. 

The  Bragg  cell  produces  several  beams,  of  which 
only  the  first  harmonic,  the  beam  that  Is  shifted  by  40 
MHz,  is  Important.  All  the  other  unwanted  beams  are  blocked 
off  by  means  of  a  beam  stop.  In  the  present  case,  a  beam 
stop  was  used  on  each  of  the  two  shifted  beams. 

Following  the  beam  stops,  beam  spacers,  TSI  Model 
9113-13,  are  used  to  for  each  set  of  beams  to  reduce  the 
distance  between  the  beam  pairs  forming  the  probe  volume. 
The  spacing  was  reduced  to  13  mm  by  using  these  prisms; 
while  maintaining  polarity,  coherence,  beam  diameter  and 
parallel  optical  paths.  The  primary  use  of  these  elements  Is 
to  control  the  fringe  spacing  in  the  probe  volume  by 
changing  the  angle  of  intersection  of  the  beam. 
Additionally,  however,  they  permit  transmission  of  the  beams 
Into  the  test  section  by  requiring  only  small  mirrors  and 
optic  components. 

The  four  beams  from  the  optic  tube  are  directed  to 
the  focusing  lens  with  five  mirrors.  The  first  two  move  the 
beams  from  10.8  cm  to  23  cm  above  the  optical  table  and  the 
next  two  reflect   the   beams  back   180°   toward   the   last 
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mirror  in  this  series.  This  mirror  directs  the  four  beams 
vertically  toward  the  transmitting  optics  lens.  The  four 
beams  then  enter  the  transmitting  optics  lens  to  focus  the 
four  beams  to  a  point.  A  final  mirror  Is  used  to  direct  this 
focused  beam  to  the  test  section.  The  focusing  lens  and  the 
last  three  mirrors  are  mounted  on  a  traversing  mechanism  to 
enable  traversing  of  the  focused  beams  in  any  of  the  three 
orthogonal  directions  as  shown  in  Figure  2.6. 
2.   Receiving  Optics  and  Photodetector 

All  the  receiving  optics  are  mounted  on  a  traversing 
mechanism  which  is  slaved  to  the  traversing  mechanism  on  the 
transmitting  side.  The  two  sides  were  aligned  to  within  0.01 
mm  and  could  be  moved  together.  A  Taskmaster  controller  is 
used  to  move  the  two  traverses  both  manually  and  through  the 
use  of  computer  programs.  See  Figure  2.5. 

The  receiving  optics  are  used  to  capture  the  light 
scattered  by  the  particles  crossing  the  probe  volume  and 
focus  this  light  onto  photodetectors.  It  consists  of  several 
components;  that  are  mounted  on  a  hardware  arm  connected  to 
the  traversing  mechanism.  The  scattered  light  is  collected 
by  a  series  of  three  lenses  that  col  1 Imate  the  light  and 
direct  It  toward  the  polarization  color  separator.  The  last 
lens  acts  as  a  spatial  filter  in  addition  to  being  a  beam 
collimator.  Its  function  Is  to  allow  light  focused  on  the 
lens  to  pass  through  while  blocking  all  the  rest  of  the 
light.  It  Is  used  to   reduce   the  stray  light  scattered  from 
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around  the  probe  volume,  as  Is  the  case  In  measurements 
close  to  a  solid  surface. 

The  polarization  beam  separator  separates  the  two 
wavelengths  present  In  the  scattered  light  into  Its  con- 
stituent colors;  green  and  blue.  This  splitter  reflects 
light  with  horizontal  polarization  upward  while  transmitting 
the  vertically  polarized  light  through  it.  This  insures  sep- 
aration of  the  two  primary  colors  in  the  scattered  light. 

Two  notch  filters  are  used  at  this  stage,  on  each  of 
the  colors  to  eliminate  residual  color  contamination.  The 
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Figure  2.5  Receiving  Optics  Schematic 
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separated  scattered  light  Is  now   focused  on  to  the  photo- 
detector  aperture  with  a  lens. 

The  photodetector  converts  the  light  signal  into  an 
electric  signal.  The  photodetector  used  is  a  photomul t lpl ler 
tube  or  PMT.  The  advantages  of  the  PMT  over  the  other  types 
are  low  noise,  high  gain  and  broad  bandwidth  while  its 
fragility  and  the  need  for  a  high  voltage  source  are  its 
disadvantages  [Ref.  2:p.  453. 

3.  Test  Section 

The  transonic  test  section  is  formed  by  inserting 
upper  and  lower  ventilated  walls  In  a  subsonic  Indraft  wind 
tunnel.  The  test  section's  dimensions  are  13  cm  by  25  cm. 


Figure  2.6  Transmitting  Side  of  Test  Section 
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A  7.5  cm  chord  NACA  0012  airfoil  was  mounted  between  the 
side  walls  at  the  centerline.  This  airfoil  Is  attached  at 
quarter  chord  to  a  external  device  to  oscillate  the  airfoil 
through  various  angles  of  attack.  The  side  walls  of  the  test 
section  are  made  of  2.5  cm  thick  plexlglas  to  allow  complete 
and  unobstructed  optical  access  to  the  airfoil.  See  Figure 
2.6.  The  mechanism  to  mount  the  model  partially  blocks  the 
views  of  the  leading  edge  of  the  model.  However,  this  was 
not  a  limitation  for  the  present  study. 
4.   Signal  Processing 

The  signal  from  the  photodetector  must  be  processed 
to  obtain  the  frequency  Information  in  the  scattered  light. 
Figure  2.7  shows  the  block  diagram  of  the  signal  processing 
components.  The  signal  from  the  PMT  is  amplified  by  an  EIN 
Model  403LA,  linear  amplifier  to  boost  signal  strength  by  37 
dB  to  transmit  the  signal  over  a  long  distance  and  to  be  of 
appropriate  strength  for  signal  processing. 

A  downmixlng  circuit  is  used  next  to  decrease  the 
effective  frequency  shift  from  40  MHz  to  around  10  Mhz  to 
Improve  the  counter  resolution.  The  two  elements  of  the 
downmixer  are  a  Tektronix  SG503  Levelled  Sine  Wave  Generator 
and  a  Hewlett-Packard,  Model  10534A,  Passive  Diode  Mixer. 
The  expected  Doppler  frequency  should  be  estimated  to  prop- 
erly set  the  mixer  frequency.  Using  the  range  of  Mach  num- 
bers expected  as  -0.2  to  +0.5,  the  free  stream  velocity 
range  can  be  calculated  using   the   set   of   equations   in 
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sect  ton  III  A  and  20°  C  as  the  ambient  total  temper- 
ature. The  velocity  range  resulting  from  this  calculation  Is 
-67  to  +167  meters  per  second.  With  the  fringe  spacing  of 
19.52  microns  and  the  second  equation  on  page  35,  the 
estimated  Doppler  frequencies  are  -3.4  to  +8.6  MHz.  Using 
the  equation  on  page  34,  the  signal  and  mixing  frequency 
should  be  expected  to  fall  within  the  range  of  36.6  to  48.6 
MHz.  By  using  a  starting  mixing  frequency  of  35  MHz,  the 
calculated  signal   frequency   reaching   the   counters  will 
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be  1.6  to  13.6  MHz.  The  mixing  frequency  should  be  adjusted 
whenever  the  tunnel  conditions  change  or  In  attempts  to 
obtain  better  counter  resolution. 

The  next  signal  processing  unit  encountered  is  the 
Series  3000  Macrodyne  LDV  Processor.  This  device  has  several 
functions:  variable  high  and  low  pass  filters,  input  and 
output  ports,  threshold  and  overload  adjustments  and  freq- 
uency counter  with  fringe  passing  selection  switch.  Both  the 
low  pass  filters,  to  remove  the  high  frequency  noise,  and 
the  high  pass  filters,  to  reduce  the  low  frequency  pedestal 
from  the  signal,  are  adjustable  to  provide  the  optimal  fil- 
tering of  the  signal.  Once  the  signal  has  been  thus  condi- 
tioned, a  frequency  counter  is  used  to  measure  the  signal 
frequency  and  convert  It  to  digital  time  data  that  can  be 
sent  to  the  data  processor.  The  counter  operates  on  two 
different  systems  to  validate  the  Doppler  signal:  a  three 
level  detection  scheme  and  a  5/8  or  10/16  comparison.  The 
three  level  detection  scheme  shown  In  Figure  2.8  must  be  met 
before  a  signal  is  validated.  A  threshold  level  is  set  on 
the  counter  to  be  Just  above  the  level  of  the  steady  state 
noise.  The  signal  must  pass  the  following  gates  to  be  ac- 
cepted: (1)  exceeds  the  threshold  setting  (2)  crosses  zero 
and  <3>  exceeds  a  negative  threshold.  To  pass  this  check, 
the  selected  number  of  fringe  crossings  must  sequence 
through  each  check  without  error.  This  detection  scheme  is 
notably   effective  in  noise  discrimination.  Once  a  signal  is 
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received  and  passes  the  first  check,  the  counter  starts  to 
measure  the  Doppler  period  by  opening  a  time  gate  and  stop- 
ping at  the  switch  set  number  of  cycles.  This  setting  of 
cycles  Is  able  to  be  selected  from  a  switch  on  the  counter 
control  panel.  Validated  data  from  the  switch  being  in  the 
10/16  cycles  position  tend  to  be  better  for  strong  Doppler 
bursts  because  the  greatest  number  of  cycles  measured  tends 
to  reduce  the  effects  of  noise  in  the  signal.  However,  this 
setting  also  restricts  the  motion  of  the  particle  through 
the  probe  volume   and   thus   Introduces  artificial  biases  of 
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the  data.  The  particle  must  traverse,  the  probe   volume   at 
its   largest   length   to  insure  validated  data  for  16  fringe 
crossings.  For  this  reason,  present  experiments  were  carried 
out  with  the  switch  set  to  the  5/8  position. 
5.   Data  Processing  Hardware 

Once  the  data  is  received  from  the  Macrodyne  signal 
processors,  it  must  be  multiplexed  and  transferred  to  the 
computer  for  velocity  calculations.  The  NASA  LDV  Computer 
Interface  is  now  used  to  control  the  multiplexing  of  data, 
sequencing  of  data,  provide  a  digital  display  to  track 
successful  counts  of  the  Doppler  frequency  and,  generally, 
control  the  flow  of  data  between  the  signal  processor  and 
the  Micro- VAX  computer.  Since  a  two  component  Laser  Doppler 
Velocimeter  Is  used  in  this  experiment,  the  multiplexer 
sends  four  channels  of  data  for  each  valid  Doppler  burst. 
The  first  data  word  contains  the  time  Information,  the 
second  contains  any  status  Information  from  the  signal  and 
the  last  two  words  contain  the  frequency  data  from  each 
component.  Due  to  the  limited  nature  of  this  experiment, 
the  first  two  words  of  data  are  not  useful  in  further 
velocity  calculations.  Therefore,  Just  the  frequencies  from 
each  particle  caused  Doppler  burst  are  used  to  calculate  the 
respective  particle  velocities. 

The  Micro-VAX  computer  stores  all  acquired  data  in  a 
buffer  set  up  by  acquisition  programs.  The  next  step  after 
data   acquisition   is   to  calculate   the   velocities.  These 
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computer  programs  reduce  the  massive  amounts  of  raw  data 
obtained  into  meaningful  velocity  data  that  can  be  output  In 
tabular  form  and  graphical  form. 

The  following  is  a  simplified  path  followed  by  the 
numerous  computer  programs  to  reduce  the  data.  At  least  1000 
samples  from  each  component  are  required  to  obtain  a  statis- 
tically steady  sample  frequency.  This  frequency  data  must 
then  be  converted  to  the  true  Doppler  frequency  by  the  fol- 
lowing equation: 

fd  =  *a  ~  fB  +  fM 

where  f^  is  the  true  frequency  resulting  from  the  particle 
traversing  the  probe  volume,  fg  is  the  signal  frequency 
measured  by  the  counter,  f3  is  the  Bragg  shift  frequency 
of  40  MHz  and  fM  is  the  frequency  of  the  downmixer.  The 
downmixer  circuit  permits  decreasing  of  the  effective  freq- 
uency shift  to  values  within  the  1  to  10  MHz  range.  The  true 
frequency  data  Is  then  averaged  over  the  samples  taken  to 
obtain  a  mean  value  and  a  sample  standard  deviation.  Using 
this  mean  value,  for  each  channel,  all  samples  that  are 
within  2  tr ,  two  standard  deviations,  of  the  mean  are  saved 
and  the  rest  are  discarded.  This  procedure  accounts  for  a- 
bout  5%  reduction  of  data  and  will  insure  samples  that  are 
too  low  or  high  will  not  bias  the  final  results.  A  new  mean 
and  standard  deviation  are  calculated  for  each  channel.  Next 
the  probe  volume   dimensions   and  parameters  are  calculated 
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with  primary  emphasis  on  the  fringe  spacing  equation: 

6  *  X  /  C2  *   sln(e/2)] 
where   *   is   the   fringe   spacing  in  microns,  X  is  the  beam 
wavelength  in  microns  and  O     Is   the  angle  of  intersection 
of  the  beams.  Knowing  the  fringe  spacing,  the   velocity   can 
be  determined  using  the  relation: 

u  =  fd  *  6 
where  u  is  the  velocity  of  the  particle  In  meters  per  second 
and   fd   is  frequency  of  the  particle  measured  by  the  count- 
ers in  MHz. 

Once  the  velocity  components  of  the  particle  are 
determined  streamwise  and  In  the  transverse  directions,  the 
computer  stores  the  results  In  a  data  table  for  further 
analysis.  The  probe  position  is  also  saved  in  this  data 
file.  The  entire  process  is  now  repeated  at  different 
locations  of  interest  in  the  flow. 

C.   LABORATORY  APPARATUS 

Besides  the  Laser  Doppler  veloclmeter  and  schl leren 
systems  described  earlier,  additional  instrumentation  was 
used  to  obtain  the  wind  tunnel  pressure  and  temperature 
data.  A  24  port  scan  1  valve  pressure  transducer  system  was 
set  up  to  measure  the  total  and  static  pressures  at  various 
points  Inside  the  operating  tunnel.  This  scanlvalve  system 
was  driven   by   a  digital  to  analog  interface,  built  on  the 
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premises,  which  allowed  computer  control  of  the  "stepping" 
and  "homing"  of  the  scanlvalve. 

Using  the  Digital  to  Analog  interface  on  the  MicroVAX 
computer,  these  signals  to  step  the  scanlvalve  system  were 
transferred  from  a  computer  program  to  the  physical  unit  for 
action.  Once  the  stepper  motor  was  sent  the  appropriate  sig- 
nal to  move,  the  Analog  to  Digital  Interface  board  on  the 
computer  sampled  the  pressure  transducer  signals  and  loaded 
the  data  Into  the  computer  for  conversion  Into  pressure  data 
and  later  analysis. 

The  atmospheric  pressure  was  recorded  from  a  mercury 
barometer  with  micrometer  scale  mounted  on  the  laboratory 
wall  near  the  wind  tunnel.  The  total  temperature  was  record- 
ed on  a  mercury  thermometer  mounted  on  the  barometer.  Both 
of  these  measurements  will  soon  be  replaced  with  Paro- 
sclentific  pressure  sensors  and  thermocouples  that  will 
have  a  direct  link  with  the  computer. 

A  TSI  Six-Jet  Atomizer,  Model  9306,  was  used  to  provide 
particle  seeding  for  the  LDV  system. 

Figure  2.9  shows  the  general  laboratory  setup  and 
location  of  the  data  processing  and  ancillary  equipment. 

D.   DATA  PROCESSING 
1  .   Introduction 

The   Laser  Doppler  velocimeter  yields  a  large  amount 
of  raw  data  during  Its   operation.   This  data  must  be  stored 
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Figure  2.9  Laboratory  Setup 
In  the  computer  In  a  logical  sequence  In  order   to  properly 
calculate  the  velocities  and  all  other  associated  data. 
2.   Computer  Programs 

The  following  discussion  pertaining  to  the  func- 
tioning of  each  of  the  computer  programs  that  was  written  to 
accomplish  the  task  of  managing  the  vast  amount  of  data  in 
the  most  efficient  and  timely  manner. 
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a.  Master  Program 

MASTER  controls  the  overall  flow  of  data  during 
each  of  the  subroutine  calls,  acquires  all  the  data  and 
reduces  it  to  a  usable  form. 

First,  the  wind  tunnel  conditions  are  sampled  by 
the  use  of  various  subroutines:  DA11J1LAB,  SETUP,  GETDATA, 
DETACH  and  OUTPUTD.  The  traversing  mechanism  is  then  posi- 
tioned with  subroutine  HWTRVT  to  move  the  probe  volume  of 
the  LDV  to  the  proper  location  in  the  flow.  Next,  the  LDV 
data  is  acquired  with  subroutine  HWWA.  The  samples  are 
taken,  velocities  calculated  using  subroutine  VELOCITY  and  a 
statistical  mean  and  standard  deviation  are  calculated.  The 
data  is  then  reduced  with  subroutine  REDUCE_DATA  to  preclude 
the  use  of  any  velocities  greater  or  less  than  two  standard 
deviations  away  from  the  mean.  Each  of  the  two  velocity  com- 
ponents are  then  recorded  along  with  the  position  in  the 
flow.  The  velocity  profiles  are  obtained  from  the  data  using 
post-acquisition  processing  programs.  An  option  is  given 
during  the  acquisition  process  to  view  a  histogram  of  the 
resulting  velocities  to  check  on  the  distribution  of  data 
samples  with  subroutine  HISTOGRAM.  A  sample  histogram  from  a 
test  run  is  shown  in  Figure  2.10. 

b.  Scan  1  valve  Programs 

First  a  signal  is  sent  to  the  scanlvalve  to 
"home"  to  port  0  for  the  first  pressure  reading.  This  is 
accomplished  by  using  the  subroutine  DA11J1LAB.  Next,  a  call 
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Figure  2.10   Sample  Histogram 

to  subroutine  SETUP  is  made  to  connect  the  computer's  in- 
ternal timing  clock  and  Analog  to  Digital  interface  board  to 
program  control.  Also,  the  sample  acquisition  rate  and  num- 
ber of  samples  are  interactively  entered. 

Now,  the  computer  is  ready  to  acquire  data  from 
the  scanlvalve  system.  Since  there  are  24  ports  to  sample  on 
the  scanlvalve,  a  loop  is  entered  in  which  samples  are 
taken,  scanlvalve  stepped  and  data  reduced  during  each  loop. 
Subroutines  GETDATA  and  DA11J1LAB  acquire  the  data  and  step 
the  scanlvalve  respectively.  DA11J1LAB  has  an  adjustable 
time  delay  to  give  time  for  the  pressure  in  each  port  to 
equalize  after  scanlvalve  stepping. 
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After  all  24  ports  are  sampled,  the  mean  and 
standard  deviation  voltages  are  calculated  for  each  port. 
The  calibration  pressure  attached  to  port  0  and  the  atmos- 
pheric pressure  at  port  1  are  now  entered  interactively  in 
consistent  units;  inches  of  mercury  in  this  case.  With  these 
two  pressures  entered,  a  calibration  line  is  calculated  and 
the  remaining  voltages  from  the  scanivalve  system  are  con- 
verted to  pressures.  Using  port  2,  settling  chamber  stagna- 
tion pressure,  and  port  7,  test  section  static  pressure,  a 
tunnel  Mach  number  Is  calculated.  See  section  III  A  for  the 
equations  used  in  these  calculations.  Next,  the  total 
temperature  is  entered  and  the  free  stream  velocity  is 
cal cu 1 ated. 

For  use  In  labelling  a  program  generated  plot- 
ting file  for  post-acquisition  analysis,  the  angle  of  attack 
of  the  test  section  airfoil  and  current  time  are  entered. 

The  last  part  of  the  scanivalve  data  acquisition 
programs  close  the  loop   on   the  analog  to  digital  Interface 
by  detaching  the  clock  and  driver.  A  diagnostic  output   file 
is  also  generated  for  post-acquisition  analysis, 
c.   LDV  Data  Acquisition  Program 

Both  the  transmitting  and  receiving  optics  of 
the  Laser  Doppler  Veloclmeter  are  mounted  on  three  axis 
motorized  traverses  that  are  slaved  together  to  move  the 
probe  volume  to  any  point  of  Interest  in  the  flow.  Using  a 
pre-existing  software   driver   for   the   Taskmaster  traverse 
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control  unit,  the  optics  could  be  moved  at   will   from   the 
computer  workstation. 

Although  a  pre-existing  software  driver  was 
available  for  the  NASA-LDV  Computer  Interface,  significant 
changes  had  to  be  made  to  the  code  for  it  to  work  on  the 
present  system.  The  main  Laser  Doppler  Velocimetry 
acquisition  subroutine  "DAMUX"  was  responsible  for  loading 
the  incoming  data  from  the  multiplexer  Into  data  arrays. 
Once  the  raw  frequency  data  Is  in  arrays,  it  has  to  be 
decoded  to  correspond  to  usable  numbers. 

The  time  word  output,  as  14  bits,  by  the  counter 
has  a  mantissa  and  exponent  format:  with  a  10  bit  mantissa 
and  a  4  bit  exponent.  The  computer  reads  it  as  a  16  bit  word 
and  adds  two  bits  that  are  not  used  for  data  handling.  These 
two  bits  are  only  used  during  direct  computer  links.  This  16 
bit  word  is  output  as  a  complementary  word  that  is  reversed; 
that  is,  the  zeros  are  output  as  ones  and  the  ones  as  zeros 
then  the  entire  word  is  reversed  (i.e.  the  least  significant 
bit  which  is  on  the  right  most  side  bit  is  moved  to  the  left 
most  bit  in  the  word).  First,  the  word  is  made  into  its  com- 
plement then  both  the  mantissa  and  the  exponent  are  reversed 
prior  to  decoding.  The  exponent  is  preset  from  the  counter 
front  panel  switch.  The  following  is  the  standard  16  bit 
word  composition  and  what   actually  results  from  the  output: 
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Correct  output. 
0111010001   0111   01    =   465x27   =  59,520 
Mantissa    Exp   NU 

Actual  output 
01   0001   0111010001 
NU   Exp    Mantissa 

Translated  by  the  computer 
0100010111  0100  01   =   279x24   =  4,464 

First  step:  Take  its  complement 
10   1110   1000101110 
NU   Exp    Mantissa 

Next:  reverse  the  order 

0111010001   0111   01   =  59,520 

Mantissa  Exp  NU 
where  NU  stands  for  Not  Used.  In  order  to  fill  the  sixteen 
bit  word,  the  14th  and  15th  bit  are  added  to  control  data 
flow  but  have  no  significance  to  the  value  of  the  data.  From 
this  last  16  bit  word  is  calculated  the  Doppler  shift  freq- 
uency. The  output  from  the  NASA-LDV  Computer  Interface  is 
actually  representative  of  the  time  elapsed  for  the  particle 
to  traverse  8  fringes  of  the  probe  volume.  The  Doppler 
period  Is  expressed  by  the  following  equation: 
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td  =  1/32  <D9D8...D0)  x  2®®   ns 

where  D9..DQ  represents  the  binary  mantissa  and  EXP  is 
the  preselected  exponent  set  on  the  counter.  From  the 
period,  the  Doppler  frequency  is  calculated  using  the 
following  equation: 

fd  =  1000/td  MHz 

This  Doppler  frequency  is  then  be  used  to  calculate  the 
velocity  of  the  traversing  particle.  Since  frequency 
shifting  is  used  to  discriminate  the  actual  direction  of  the 
particle  through  the  probe  volume,  the  Doppler  frequency 
must  be  adjusted  to  account  for  the  Bragg  and  mixing 
frequencies.  See  equation  on  page  34  for  this  relationship. 
Once  the  actual  frequency  of  the  particle  is  determined,  the 
fringe  spacing  is  calculated  by  the  first  equation  on  page 
35  and  the  velocity  is  found  using  the  second  equation  on 
the  same  page. 
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d.   LDV  Support  Programs. 

Three  subroutines  are  used  to  provide  the 
supporting  functions  to  the  LDV  data  acquisition.  These  are 
VELOCITY,  HISTOGRAM  and  REDUCE_DATA. 

VELOCITY  is  responsible  for  calculating  the 
fringe  spacing,  extracting  the  particle  frequency  from  the 
signal  and  then  calculating  the  two  velocities. 

HISTOGRAM  takes  the  velocity  data  and  plots 
a  histogram  of  the  percentage  of  velocities  that  fall  within 
a  given  range.  The  x  axis  of  the  histogram  Is  divided 
equally  Into  8  segments  from  the  minimum  to  maximum  velocity 
values.  The  option  to  plot  a  velocity  histogram  is  given 
after  each  LDV  acquisition  to  spot  check  the  distribution  of 
the  data.  A  nearly  Gaussian  or  normal  distribution  is 
expected  if  the  system  is  operating  properly.  This  will  also 
show  the  effects  of  noise  as  evidenced  by  the  broad  band  of 
the  histogram  values. 

REDUCE_DATA  reduces  the  velocity  data  to 
only  those  values  falling  within  two  standard  deviations  of 
the  mean.  This  results  in  a  4.5%  rejection  of  data  which 
will  reduce  the  effect  of  high  frequency  noise  on  the 
counters  and  tend  to  reduce  the  effect  of  very  large 
particles  moving  slowly  through  the  probe  volume. 
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III.  DATA  ACQUISITION 

A.   TUNNEL  CONDITIONS 

In  order  to  obtain  good  correlation  between  the  schlle- 
ren  photographs  and  the  Laser  Doppler  Veloclmeter  obtained 
data,  the  operating  tunnel  conditions  would  have  to  be 
nearly  Identical.  Two  sets  of  conditions  were  chosen:  one 
corresponding  to  an  approximate  Mach  number  of  0.3  and  the 
other  0.5.  Once  the  tunnel  was  running  at  operating 
conditions,  the  computer  program  MASTER  calculates  the  test 
section  Mach  number  and  free  stream  velocity  according  to 
the  following  relationships.  Assuming  Isentroplc  flow  In  the 
test  section,  the  section  Mach  number,  M,  Is  calculated  from 
the  stagnation  pressure,  PQ,  and  the  static  pressure,  P, 
using  the  following  relationship: 

M  »  <([P0/P](lr-1/r)-l)x<2/[y-l]>)1/2 

where  y  =  1.4.  The  resulting  calculated  Mach  numbers  for 
the  two  conditions  are  as  follows:  0.295  and  0.450.  From 
these  Mach  numbers,  the  corresponding  free  stream  velocities 
are  calculated. 

V  -  M  *  a 

where  a  Is  the  sonic  velocity. 

a  -  <yRT>1/2 

where  R  Is  the  universal  gas  constant  and  T   Is   the   static 
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temperature  In  the  test  section.  The  static  temperature  can 
be  calculated  from  the  total  temperature,  TQ,  measured  on 
the  wall  near  the  tunnel  by  the  next  relationship: 

T  =  Tq/U  +  [(y-l)/23*M2) 

B.   SCHLIEREN  PHOTOGRAPHS 

1 .  Introduction 

Once  the  wind  tunnel  was  running  at  the  first  of 
two  operating  conditions,  manometer  board  readings  were 
recorded  and  Mach  number  calculated  using  the  first 
equation  on  the  preious  page.  The  schlleren  photographs  were 
then  taken  at  four  different  angles  of  attack  and  three 
different  strobe  intensities  to  obtain  the  optimal 
photographs. 

2.  Test  Conditions  and  Procedure 

Using  the  procedures  detailed  in  section  II  A,  the 
schlleren  photographs  were  taken  at  the  settings  and 
conditions  listed  in  Table  I,  below. 

Upon  completion  of  the  wind  tunnel  runs,  the  35mm 
film  was  developed  and  printed.  From  these  schlleren 
photographs,  the  conditions  resulting  In  the  best  depiction 
of  the  wake  profiles  were  chosen  to  represent  each  Mach 
number  and  angle  of  attack.  These  are  annotated  with  an 
asterisk,  *,  In  Table  I,  below. 

Using  these  photographs,  a  short  HP-41  CX  calculator 
program  was  written  to   relate   photographic   dimensions   to 
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TABLE  I.   SCHLIEREN  PHOTOGRAPH  PARAMETERS 


Strobe 

a    (deg) 

Intensity  (Hz) 

0.0 

20 

0.0 

70 

0.0 

150 

7.7 

20 

7.7 

70 

7.7 

150 

15.5 

20 

15.5 

70 

15.5 

150 

23.6 

20 

23.6 

70 

23.6 

150 

0.0 

20 

0.0 

70 

0.0 

150 

7.7 

20 

7.7 

70 

7.7 

150 

15.5 

20 

15.5 

70 

15.5 

150 

£1 

#  0.30 
0.30 
0.30 
0.30 
0.30 

#  0 .30 

*  0 .30 
0.30 
0.30 
0.30 
0.30 
0.30 

*  0.45 
0.45 
0.45 

*  0.45 
0.45 
0.45 
0.45 

#  0.45 
0.45 


real  test  section  dimensions.  Once  the  real  dimensions  of 
the  wake  were  known,  the  Taskmaster  traverse  mechanism  was 
calibrated.  This  yielded  a  1:1  relationship  in  the  chordwise 
direction  and  a  1:4  relationship  in  the  vertical  direction. 
This  means  that  the  Taskmaster  digital  control  panel  reads 
true  values  for  the  X  direction,  chordwise,  and  four  times 
the  real  value  in  the  Z  direction,  vertical.  Arbitrarily 
choosing  the  trailing  edge  of  the  airfoil  at  0°  «  as  an 
Initialization  point  for  the  Taskmaster,  a  table  of  LDV 
sampling  points  was  made  for  the  first  three  angles  of 
attack.  See  Table  II,  below,  for  a  sample  of  the  Taskmaster 
table  at  a  -  7.7° 
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Table  II.   TASKMASTER  LDV  SAMPLING  LOCATIONS 


Points 

I 

--> 

Alpha 

=  7.7  deg 

20  2.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  - 

-1.5 

-2 

-2.5 

20  1.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  - 

-1.5 

-2 

-2.5 

20  1.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  - 

-1.5 

-2 

-2.5 

20   .50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  - 

-1.5 

-2 

-2.5 

20   .00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5  - 

-1.5 

-2 

-2.5 

19  -.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

0 

-.5 

-1.5 

-2 

-2.5 

17  -1.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

-1.5 

-2 

-2.5 

17  -1.50 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

1.5 

1 

.5 

-1.5 

-2 

-2.5 

11  -2.25 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

-2.5 

12  -3.00 

7 

6.5 

6 

5.5 

5 

4.5 

4 

3.5 

3 

2.5 

2 

-2.5 

176   (Total 

locations  saapled) 

C.   LASER  DOPPLER  VELOCIMETRY  DATA  ACQUISITION 
1 .   Calibration  of  the  LDV  System 

The  relationship  between  the  Doppler  frequency  and 
the  velocity  of  the'partlcle  traversing  the  probe  volume  Is 
given  by  the  second  equation  on  page  35.  The  particle 
velocity  Is  directly  proportional  to  the  Doppler  frequency 
by  a  factor  of  the  fringe  spacing.  By  the  first  equation  on 
page  35,  fringe  spacing  la  only  a  function  of  the 
transmitting  lens  angle  and  wavelength  of  the  light.  Since 
the  wavelength  of  the  light  Is  accurately   known   to  within 


48 


0.05  microns,  or  0.01%,  an  accurate  lens  angle  Is  all  that 
is  needed  for  calibration  [Ref.  43. 

The  following  relationship  was  used  to  calculate 
the  lens  angle: 

e  =  2  *  tan'1  (E  *  d  /  (2  *  f)) 

where  E  is  the  beam  expansion  ratio,  1.0  in  this  case:  d  is 
the  incident  beam  spacing,  preset  at  13  mm  with  TSI  Model 
9113-13;  and  f  is  the  wavelength  of  the  light. 

The  lens  angle  in  this  case  was  not  accurate 
because  the  focal  length  of  the  focusing  lens  was  not 
known.  Therefore,  calibration  was  necessary  to  Insure 
accurate  velocity  measurements  were  obtained. 

The  calibration  method  chosen  was  to  operate  the 
wind  tunnel  at  known  conditions,  calculate  the  free  stream 
velocity  from  the  scanivalve  pressure  data  and  then  move 
the  probe  volume  via  the  Taskmaster  to  coincide  with  a 
position  in  which  the  free  stream  velocity  could  be 
measured.  From  this  procedure,  a  constant  could  be  applied 
to  the  velocity  equation  as  follows: 

u  =  fd  *  6  *  (Calibration  Constant) 

This  calibration  constant  will  compensate  for  the  estimated 
focusing  lens  focal  length. 

2.   Test  Conditions  and  Procedure 

Once  calibration  is  completed,  the  tunnel  is 
brought  to  the  test  conditions,  the  Taskmaster  is  moved   to 
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the  first  sample  point  and  the  MASTER  computer  program  Is 
run  with  the  adjusted  fringe  spacing  to  reflect  the 
cal lbrat ion . 

At  each  sampling  point,  1000  samples  will  be  taken 
to  Insure  statistical  accuracy  of  the  mean  flow  calcula- 
tions. The  scheduled  sampling  points  can  be  seen  In  Table 
II,  page  48.  After  each  Taskmaster  traverse,  the  oscillo- 
scope Is  checked  to  insure  no  drastic  loss  of  the  signal 
occurred  during  the  movement.  If  the  Doppler  signal  de- 
creases, slight  adjustments  to  the  receiving  optics  traverse 
may  be  needed  to  regain  the  signal  strength.  To  optimize  the 
data  received,  slight  adjustments  may  have  to  be  made,  also, 
to  the  selected  mixing  frequency,  mixing  frequency  ampli- 
tude, counter  threshold  setting  or  high  and  low  pass  filter 
sett  1 ngs. 

After  Table  II  data  points  are  sampled,  the  airfoil 
angle  of  attack  or  test  Mach  number  is  changed  and  the  above 
procedure  is  repeated. 
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IV.     RESULTS    AND    DISCUSSION 

A.      SCHLIEREN 

The  first  thesis  goal  was  attained  early  enough  to 
provide  qualitative  information  concerning  the  airflow 
around  and  in  the  wake  of  the  airfoil.  Schlleren  35mm 
photographs  were  produced  for  both  selected  mach  numbers, 
0.3  and  0.45,  and  for  four  different  angles  of  attack;  0°, 
7.7°,  15.5°  and  23.6°.  The  best  photographs  for  each 
Mach  number  and  the  first  three  angles  of  attack  are  dis- 
played In  Appendix  A,  pages  55-57.  From  these  photographs, 
LDV  sampling  points  were  chosen  as  seen  from  Table  II  on 
page  48. 

Three  turbulent  areas  are  evident  on  all  the  schlleren 
photographs.  As  seen  in  Figure  A. 4  on  page  56,  the  airflow 
Is  slightly  disturbed  upstream  of  the  airfoil  due  to  the 
presence  of  a  pitot-static  tube  extending  into  the  airflow 
Just  forward  of  the  airfoil.  This  tube  has  been  subsequently 
removed  to  clean  up  the  flow.  At  the  top  of  the  same  photo- 
graph, a  wedge  shaped  disturbance  can  be  seen.  This  and  the 
lower  floor  turbulent  layer  are  results  of  the  transonic 
wall  slottings  found  on  both  the  roof  and  the  floor  of  the 
test  section.  Although  a  concern  for  later  unsteady  flow 
studies,  these  two  areas  do  not  Interfere  with  the  wake 
profiles  resulting  from  a  steady  flow  analysis. 
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B.   LASER  DOPPLER  VELOCIMETER 

Due  to  equipment  problems  and  time  constraints,  the  LDV 
data  acquisition  part  of  this  thesis  was  not  able  to  be 
accomplished.  All  systems  are  now  operational,  to  include 
computer  programs,  however,  there  is  inadequate  time  to 
acquire  data  and  analyze  the  velocity  profiles  obtained. 

Calibration  runs  were  made  with  a  limited  degree  of 
success  when  the  wind  tunnel  was  run  under  Mach  0.3  and  the 
Laser  was  not  powered  over  2  watts.  The  problem  with  the 
Laser  is  easily  repaired,  however,  time  consuming.  The 
upper  limit  on  the  wind  tunnel  Mach  number  is  caused  by  the 
evaporation  of  the  seeding  particles  in  high  speed  flow. 
This,  too,  can  be  easily  overcome  by  the  proper  selection 
of  seeding  particles  and  related  atomizer. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  thesis  provided  Invaluable  assistance  to  the  dynamic 
stall  project  at  the  Fluid  Mechanics  Laboratory.  The  schlie- 
ren  photography  was  able  to  provide  the  qualitative  data  on 
the  steady  state  flow  field.  With  this  accomplished,  the 
schl ieren  system  in  this  laboratory  was  able  to  be  easily 
tailored  to  the  next  step:  unsteady  flow  analysis. 

The  LDV  work  provided  the  laboratory  with  the  essential 
computer  programs,  hardware  Interfaces  and  drivers,  optics 
systems  and  auxiliary  required  equipment  to  enable  the  pro- 
ject to  proceed  with  very  little  additional  work  to  achieve 
the  steady  state  velocity  profiles.  The  ground  work  has  been 
accomplished  to  enable  the  project  to  continue  forward  with 
the  systems  on  hand  to  study  the  effects  of  dynamic  stall  in 
both  the  steady  and  unsteady  environments. 

B.  RECOMMENDATIONS 

The  schl ieren  system  is  fully  functional,  however,  the 
Laser  Doppler  Velocimeter  system  still  requires  a  little 
additional  effort  to  make  it  productive.  The  Laser  doughnut 
condition  needs  to  be  addressed  to  enable  higher  powers  to 
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be  utilized.  The  types  of  seeding  ayailable  to  the  system 
needs  to  be  addressed  to  insure  operabillty  of  the  tunnel 
and  data  acquisition  systems  above  Mach  0.3.  With  these 
minor  fixes  and  a  little  fine  adjustments  of  the  signal 
processing  equipment,  the  LDV  is  ready  for  action. 
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APPENDIX  A:  SCHLIEREN  PHOTOGRAPHS 


Figure  A.l   Mach  Number  =0.3,  a   =  0' 


Figure  A.  2  Mach  Number  =  0.45,  «  =  0° 
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Figure  A. 3  Mach  Number  =  0.3,  a  «  7.7° 


Figure  A. 4  Mach  Number  =  0.45,  a   ■  7.7° 
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Figure  A. 5  Mach  Number  =  0.3,  a   =  15.5° 


Figure  A.  6  Mach  Number  =  0.45,  a   =  15.5° 
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APPENDIX  B:  COMPUTER  PROGRAMS 


PROGRAM  MASTER 


C 

C  Master  Program  written  by  Maj  Richard  R  Ryles.  NPS  Student  28  Apr  87 

C  Last  Update  6  Sap  87 

C 

C  This  program  will  gather  dota  fro*  the  scanivalve  system,  move  the 

C  Taskmaster  and  acquire  LOV  velocity  data. 

C   This  progrom  should  be  used  with  subroutine  SETUP,  GETDATA.  DETACH,  OUTPUTD 

C   HWTRVT.  HWWA.  HISTOGRAM  and  VELOCITY. 

C — 

INCLUOE    ,SYS$LIBRARY:LIOSET.FOR• 


C 

C       Oeclare  variables 
C— 


INTEGER  DEV_ID.CLK_ID 

INTEGER  NCHANL.HISTO.BLUE.GREEN 

INTEGER  NPT.NPTS.NSAMPLES. COUNT 

INTEGER    I. J. NSTEPS. I ER. ANSWER. END. TALLY, TIME 

REAL  VOLTAGES(5eee).U_U(2ee).V_V(2ee).U_V(2ee).DELX.DELZ 

REAL'S      SUM.AMEAN.X.XX.Y.YY.STD_OEV.STO_DEV_V.STD_DEV_U.SUM_V,SUM_U 

REAL'S      SUM  VERT,SUM_HORIZ 

REAL  RATE.SCANI_VXX.T(e:23.2).SCANI_PRESS(e:23).VERT_VEL(82ee) 

REAL  CALIB  PRESS. SLOPE.  INTERCEPT. HORIZ_VEL(829e)  .  FINAL_V.FINAL_U 

REAL  X   POSITION. Z_POSITION.MEAN(2).M.FREE_STREAM. TO. ATMOS  PRESS 

real       x_LOC(2ee)  .z_LOC(2ee)  ,u(2ee)  .v(2ee)  .alpha,  zero,  corr.factor 

REAL'S     UU.W.UV 
CHARACTER* 1    ANS.VOLT_ANS 
LOGICAL  HOME 


C 

C       Scanivalve  dota  acquisition 
C- 


OP EN ( UN I T-99 . F I LE- ' RESULTS . OAT *. STATUS- ' NEW' ) 

TALLY  -  S 

HOME  -  .TRUE. 

NSTEPS  -  e 

CALL  DA11J1LAB  ( HOME. NSTEPS. I ER) 

IF  (IER.NE.1)  GOTO  17» 

ANSWER  -  6 

WRITE  (•••)  'Skip  to  Taskmaster?  (y/n)  [DEF  -  No]* 

READ  (».1)  ANS 

IF  (ANS.EQ. ,Y'.OR.ANS.EQ.,y')  GOTO  9t 

WRITE  (•••)  'Print  the)  scanivalve  diagnostic  data  to  o  file?  (Y/N)'. 
ft  '  (DEF  -  No]' 

READ  (».1)  VOLT_ANS 

CALL   SETUP    (CLX_ID.DEV_ID.RATE.NCHANL.NPT.NPTS) 
16  CONTINUE 

HOME  -      FALSE. 

NSTEPS  -    1 

DO   56   J   -  6.    23 
WRITE    (».2)    J 

IF    (J.NE.e)   CALL  DA11J1LAB    ( HOME. NSTEPS. I ER) 
IF    (IER.NE.1)    GOTO   15s) 
20  CONTINUE 

CALL  GETDATA    (DEV.ID.NPT .NPTS.  VOLTAGES) 

C    •    •    • eeeseeeeeeeee 

C  Calculate   neon   and   standard   deviation 


C 


sum  -  e  e 

DO  30    I    -   2.    NPT 

SUM  -  SUM  +  VOLTAGES(I) 
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36       CONTINUE 

AMEAN  -  SUM  /  FLOAT (NPT-1) 

SUM  -  8. 08 

DO  48  I  -  2.  NPT 

X  -  VOLTAGES(I)    -  AMEAN 
Y  -  X   •   X   •    (1    8/FL0AT( NPT-1)) 
46  SUM  -  Y  ♦   SUM 

STD.DEV  -  SORT (SUM) 

IF    (VOLT.ANS.EQ.'y'.OR.VOLT.ANS.EQ.'Y')   CALL  OUTPUTD    (RATE. NPT. 
k  NCHANL.  AMEAN. STD.DEV, J, VOLTAGES) 

SCAN I .VOLT f J.I)    -  AMEAN 
58  SCAN I .VOLT (J. 2)   -  STD  DEV 

CALL  DETACH   (CLK_ID.DEV_ID) 

IF(VOLT_ANS.EQ. 'Y' .OR. VOLT.ANS.  EQ. 'y')    THEN 

OPEN (  UN  I T- 1 5 . F I LE- ' SCAN I.VOLT . DAT ' . STATUS- ' NEW ' ) 
DO  68  I  -  6.   23 
66  WRITE   (15.3)    I.    SCANI_VOLT(I . 1 ) .    SCAN I _VOLT ( I . 2 ) 

ENOIF 

C •••••••••••••• 

C            Calculate   scanivalve   calibration    line 
C            and   convert   all    voltages   to  pressure 
C   •   •   •   e   e   •   e ••••••••• 

WRITE   (•••)    'Enter  Scanivalve  Calibration  Pressure:' 

READ   (•.•)   CALIB.PRESS 

WRITE(» .•)  'Enter  Calibration  Pressure  Correction   Factor:' 

READ( • . • )CORR_FACTOR 

CALIB.PRESS  -  CALIB.PRESS  +  CORR_ FACTOR 

WRITE(« .•) 'Enter  atmospheric  pressure:' 

READ(«.»)   ATMOS.PRESS 

SLOPE  -    (ATMOS_PRESS-CALIB_PRESS)/(SCANI_VX)LT(1l1)-SCANI_VOLT(e.1)) 

INTERCEPT  -  CALIB.PRESS  -  SCANI_VOLT(e. 1 )»SLOPE 

WRITE(99.«) 'CALIB.PRESS'. CALIB.PRESS 

WRITE(99.»)'CORR_FACTOR'  .CORR.FACTOR 

WRITE(99. •) 'ATMOS.PRESS' .ATMOS.PRESS 

OPEN  (UN  IT- 16.  FILE- 'SCAN  I  PORT.  DAT '.STATUS- 'NEW') 

DO  76  I  -  6.  23 

SCANI.PRESS   (I)   -  SCANI.VOLT(I.I)    •   SLOPE  +   INTERCEPT 
78  WRITE(16.«)I.    SCANI.PRESS(I) 

CLOSE(16) 

IF   (VOLT.ANS.EQ.'Y'  .OR.  VOLTJYC.  EQ.  > '  )    THEN 
DO  86   I   -  6.   23 
88  WRITE   (15.4)    I. SCAN I .PRESS (I ) .SCANI_VOLT(I . 2)»SLOPE+INTERCEPT 

END  IF 

CL0SEO5) 

WRITE(«.«)'    ' 

WRITE(«.«) 'Equation  of    the  calibration    line  -' 

WRITEf».»)'Y  -'.slope.'    X  ♦  *. intercept 

WRITE(99.«)'Y  -'.slope.'    X  ♦'.intercept 

WRITE(«.»)'    ' 

C   .......... ••••• 

C  Calculate    tunnel   Moch   and   velocity 

C •••• 

IF(SCANI.PRESSm.EQ.e.e)    SCANI    PRESS(7)   -    1.8 

IF(SCAN1_PRESS(2).LT.SCANI_PR£SS(7))   GOTO  85 

M  -  SORT(((SCANI_PRESS(2)/SCANI_PRESS(7)).«(2/7.)-1.)«5.) 
85  WRITE(*.*)'Tunnel   Mach  Number  -'.M 

WRITE(99.«) 'Tunnel   Mach  Number  -' ,M 

WRITE(».»)'    ' 

WRITE(«.*)'Enter  Total    Temperature    (C)' 

READ(«.«)    TO 

WRITE(99.») 'TOTAL   TEMP'.    TO 

WRITE(«.*)'    ' 

TO  -   TO  ♦  273.15 

FREE.STREAM  -  M  •   SORT ( 481    52«T0/(1 .+(M«»2/5. ))) 

WRITE(«.«)'Free  Streom  Velocity   (mps)   -'.FREE.STREAM 

WRITE(99. «)'Free  Streom  Velocity    (mps)   -'.FREE.STREAM 

I F( FREE.STREAM. EG. 8. 8)    FREE.STREAM  -    1.E-18 

WRITE(».»)'    ' 
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WRITE(e.e)'Entor  ongle  of  ottoek  (dog):' 

REAO(e.e)  ALPHA 

WRITE(99.e)' ALPHA".  ALPHA 

WRITE(« .•) 'Enter  Current  Tilts:  (4  digit  intsgsr)' 

Rsod(e.e)  TIME 

WRITE(99.eWlME'.TIME 

WRITE  (•.•)  'Skip  to  end?  (y/n)  [OEF  -  No]' 

READ  (».1)  ANS 

IF  (ANS.EO.  'y'  .OR.ANS.EQ.  'Y')  GOTO  169 

C 

C       Taskmostsr  mover  and  LOV  data  acquisition 

C 

96    WRITE  (•.•)  'Which  color  it  connected  to  ths  loft  (#3)  channol?'. 
k  •    (Bluo-I.Grssn-9)' 

READ  (*,9)  BLUE.GREEN 

WRITE(99.e)'BLUE_GPEEN' ,BLUE_GREEN 

IF(FREE_STREAM.EQ.9.9)    FREE_STREAM  -   1  .  E-6 

OPEN  (UNIT-25. FILE-" FINAL. OAT' , STATUS-' NEW' ) 

WRITE  (25. •)    •  Taskmaster  (Final. Oat)' 

WRITE  (25. •)  '      pos  (in)     vslocity/std  dov   vslocity/std  dov' 

WRITE  (25.16) 

WRITE  (•.•)  'Skip  to  LDV?  (y/n)  [Dof  -  NO]' 

READ  (».1)  ANS 

IF  (ANS. EQ. 'Y' .OR.ANS.EQ. 'y')  GOTO  119 

IF(TALLY.EO.e)  THEN 

WRITE(«.«) 'Entsr  X-Position  for  calibration' 
READ(«.«)  X.POSITION 

WRITE(». •) 'Entsr  Z-Potition  for  calibration' 
READ(».»)  Z.POSITION 

ENDIF 
196   CALL  HWTRVT  (DELX.DELZ) 

WRITE  (•.•)  'Movs  the  Taskmaster  again  prior  to  LDV  Acquisition?  (Y/N)' 

READ  («.1)  ANS 
Ceeeeeeeeeeeeeeeeeeeeee 

C     Record  position  of  LDV  at  all  times 

C  •  •  •  •  •  • ••••••••••• 

IF(ANS.EQ. 'Y*. OR.ANS.EQ. *y')  THEN 
X  POSITION  -  X_POSITION  +  DELX 

Z  POSITION  -  Z.POSITION  +  OELZ 

goto  iee 

ENDIF 

C  •  •  • •••••••••• 

C  Acqui  re    LDV  dato 

Ceeeeeeeeeeeeeeeeeeeeee 

1 10        CALL   HWWA    (NSAMPLES. ANSWER .BLUE.GREEN) 
C    •    • •••••••••••• 

C      Increment  position  and  record  it 

C     Read  velocity  data,  sun  it  and  find 

C        standard  deviation 

C  •  •  •  •  •  • • ••• 

X  POSITION  -  X.POSITION  +  DELX 
Z  POSITION  -  Z_POSITION  +  OELZ 
SUM.VERT   -  9.D9 

sum.horiz  -  eoe 

OPEN    (UNIT-35.FILE--VEL0CITY.DAT- . STATUS- ' OLD ' ) 

DO    129    I    -   1.    NSAMPLES 

IF    (BLUE  GREEN.  EQ.9)    READM5.7)    VERT_VEL(  I  )  .HORIZ.VELfl  ) 
IF    (BLUE.GREEN. EQ.l)    REA0(35.7)    H0RIZ_VEL( I ) . VERT_VEL( I ) 
SUM.VERT  -   SUM  VERT   +   VERT.VEL(I) 
129  SUM.HORIZ  -   SUM.HORIZ   ♦   HORIZ.VEL(I) 

CL0SE(35) 

FINAL   V  -   SUM.VERT   /   FLOAT (NSAMPLES) 

FINAL.U   -   SUM.HORIZ   /    FLOAT (NSAMPLES) 

SUM.V   -   0.D9 

SUM.U   -   6.M 

uu  -  e  oe 
w  -  e  09 
uv  -  e  09 
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00  1M  I  -1,  NSAMPLES 

X  -  VERT_VEL(I)  -  FINAL.V 
XX  -  H0RIZ_VEL(I)  -  FINALJJ 
Y  -  X  •  X  •  (1.6/FLOAT(NSAMPLES-1)) 
YY  -  XX  •  XX  •  (1.e/FL0AT(NSAMPLES-1)) 
UU  -  UU  +  ((HORIZ_VEL(I)"2)-(FINALJ>«2)) 
W  -  W  ♦  ((VERT_VEL(I)««2)-(FINAL_V«*2)) 
UV  -  UV  +  ((HORIZ_VEL(I)«VERT_VEL(I))-(FINAL_U»FINAL_V)) 
SUM_V  -  Y  +  SUM.V 
13B      SUM_U  -  YY  +  SUM_U 

UU  -  UU  /  FLOAT (NSAMPLES) 

W  -  W  /  FLOAT(NSAMPLES) 

UV  -  UV  /  FLOAT (NSAMPLES) 

STO  DEV_V  -  SORT  (SUM.V) 

STD_0EV_U  -  SORT  (SUMJJ) 

WRITE(«.«)  '  ' 

WRITE  (•.•rveloci  t  let  and  std  devs  -' 

WRITE   (•.•)'U  di r'.FINAL_U,STD_DEV_U 

WRITE   (•••),V  dir\FINAL_V.STD_DEV_V 

WRITE   (99.»)'Velocit ies   and   std  devs  ■' 

WRITE   (99.»)'U  dir'.FINAL_U.STD_DEV_U 

WRITE   (99.0'V  dir' ,FINAL_V.STD_DEV_V 

WRITE   (•••)•     ' 

WRITER.  O'Fi  rst    cut   RMS   values' 

WRITE(« . •  ) 'U-RMS.V-RMS' .UU. W 

WRITE(«.»)'UV-RMS'.UV 

WRITE(99.»)'First   cut   RMS   values' 

WRITE(99.*) 'U-RMS.V-RMS' .UU.W 

WRITE(99.«)'UV-RMS'.UV 

C  •  •  •  •  •  • • •••• 

C     Eliminate  data  outside  two  sigmo  from  mean 
C  •  •  •  • •••••••••••• 


IF  (BLU€_GREEN.E0.1)  THEN 

MEAN(1)  -  FINALJJ 

MEAN(2)  -  FINAL_V 

CALL  REDUCEJ5ATA  (MSAMPLES. COUNT. END. MEAN. STD.DEVJJ. STD  DEV.V. 
!  BLUE.GREEN) 

FINAL_U  -  MEAN(1) 

FINAL_V  -  MEAN(2) 

WRITE  (99. •) 'Reduced  velocities  and  std  devs  -' 

WRITE  (99.»)'U  dir' ,FINAL_U.STD_DEV_U 

WRITE  (99.«)'V  dir' .FINAL.V.STD  DEV.V 


ELSE 


ENOIF 


MEAN(2)  -  FINALJJ 
MEAN(1)  -  FINAL  V 
CALL  REDUCE.DATA  ( NSAMPLES. COUNT. END. MEAN. STD_DEV_V.STD_DEV_U. 

BLUE.GREEN) 
FINALJJ  -  MEAN(2) 
FINAL  J/  -  MEAN(1) 

WRITE  (99. •) 'Reduced  velocities  and  std  devs  ■  ' 
WRITE  (99.«)'U  dir'.FINAL_U.STD_DEV_U 
WRITE  (99.e)'V  di r ' . FINAL.V. STD_DEV_V 


C     Calculate  RMS  values 

C  ..•*.••*«*.•..  . 

OPEN(UNIT«43 . FI LE-'VELOCI TY . DAT ' . STATUS- 'OLD ' ) 
UU  -   8   D0 

w  -  e.oe 
uv  -  e.oe 

DO  216   I    -   1.    END 

IF   (BLUE_GREEN.EO.e)    READ (43. 7)    VERT   VEL( I ) .HORIZ_VEL( I ) 
IF  (BLUE.GREEN.EO.I)   READ(43,7)    HORIZ_VEL( I ) ,VERT_VEL( I ) 
UU  -  UU  +    ((HORIZ_VEL(I)»«2)-(FINAL_U««2)) 
W  -  W  +   ((VERT_VEL(I)«»2)-(FINAL_V.«2)) 
210  UV  -  UV  +   ((H0R1Z_VEL(I)»VERT_VEL(I))-(FINAL_U»FINAL_V)) 

CL0SE(43) 

UU  -  UU  /  FLOAT(END) 
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W  -  W  /  FLOATED©) 
UV  -  UV  /  FLOAT(END) 
WRITE(«.«V    ' 

WRITE(«.») 'Reduced  RMS  values' 
WRITEf-.O'U-RMS.V-RMS'.UU.VV 
WRITE!  •.•)'UV-RMS\UV 
WRITE(99.«)'U-RMS.V-RMS' .UU.W 
WRITE(99.»rUV-RMS'  .UV 

WRITE   (25.8)   X_POSITI0N.Z_POSITI0N.FINAL_U.STD_DEV_U.nNAU_V.STD  DEV  V. 
k  UU.W.UV 

TALLY  -   TALLY  +    1 

C   • • •••• 

C            Options    to  view  histogram  and   toke  dato  again 
C seeeeeeeeeeeeeeee 

WRITE(«.»)'     ' 

WRITE  (•••)  'Do  you  wont  a  Histogram?  (1/6)' 

READ  (•.•)  HISTO 

IF  (HISTO. E0.1)  CALL  HISTOGRAM (COUNT. END. BLUE  GREEN) 

WRITE  («.6) 

READ  (».1)  ANS 

IF  (ANS.EO. 'n'.OR.ANS.EO. 'N')  GOTO  146 

GOTO  iee 
146   CLOSE  (25) 

CL0SE(99) 
c  ............  • .  .  .  . 

C     Convert  output  data  into  proper  format  for 

C     "PLOTCV"  Vector  Plot  of  velocity  profiles 

C ....••••••••••••••• 

OP EN (UN I T-45. FILE- 'FINAL. DAT '.STATUS- 'OLD') 

READ(45.14) 

OP EN ( UN I T-46. FILE- 'GRA. DAT '.STATUS- 'NEW') 

WRITE(46.13)  M.FREE.STREAM. ALPHA, TIME. 11.1 6. It 

ZERO  -  6.6 

DO  15t  1-1. TALLY 

READ(45 . 1 1  )X_LOC( I ) . Z_LOC( I ) ,U( I ) ,V( I ) ,U_U( I ) . V_V( I ) .U_V( I ) 

X_L0C(I)-X_L0C(I)+3.t 

Z_LCC( I )-(Z_LOC( I )+2 . 5)« . 24 

U(I)  -  U(I)  /  FREE.STREAM 

V(I)   -  V(I)    /   FREE.STREAM 

U_U(I)   -  U_um   /  FREE.STREAM 

V_V(l)  -  V_vh)  /  FREE_STREAM 

U_V(I)   -  U_V(I)   /   FREE_STREAM 


15t  WR 

k 
CL0SE(45) 
CLOSE (46) 

166        WRITE    (».5) 


ITE(46. 12) I ,X_LOC( I ) .Z_LCC(I ) .U( I ) ,V( I ) .ZERO.U_U( I) ,V_V( I ) .ZERO. 
U_V( I ) . ZERO 


C  Format    statements 

C- 


1  FORMAT  (A1) 

2  FORMAT  (1X. 'Sconivolve  f    '.12) 

3  FORMAT  (iX.'Port  #'. 12. 3X. 'Average  Voltage  -'.F12.4.3X. 

*  'Standard  Deviation  -'.F12.4) 

4  FORMAT  (1X.  Port  #'. 12. 3X. 'Pressure  -' .F12. 4, *  (psi)'.3X. 
k  'Standard  Deviation  -' .F12. 4. '  (psi)') 

5  FORMAT  (//IX.'This  is  the  ENDl '//) 

6  FORMAT  (/IX.Move  taskmaster  again  and  take  data?  (y/n)  [DEF  -  Y]'/) 

7  FORMAT  (1X.2F11  .4) 

8  FORMAT  (1X.2F7.2.F16 .  2.FB.2.F16. 2.F8.2. 1P3E12.4) 

9  FORMAT  (II) 

11  FORMAT (1X.2F7  2. F16 . 2  ,8X . F16 . 2 ,8X . 1P3E12 . 4) 

12  FORMAT(1X.I3.1P18E12.4) 

13  FORMAT (IX.    GRA. DAT.    MACH   NO   -    '.F4.2.'.    FREE    STREAM  VELOCITY   -    '.F5.1. 

*  '.    ALPHA  -    ' .F4. 1 ,/1X. 'September   9.    1987' , 4». i4. '    hrs.' 

ft  //4X.NPR0F      NSAMPL        NCOEF '/7X . 12 . 6X . 12 . 6X . I2/1X . '    ID  X-POS'. 

k  Z-POS  U  V  W  Ut2'. 

*  Vt2  Wt2  UV  UW') 
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14    FORMAT  (//) 

16    FORMAT(6X.•x^6X.,2,.12X,•u•,17X.•v^11X.•ut2,.9X,•vt2,.1eX.,uv•) 
STOP 

C 

C       Sconivolv*  Error  Stotoiionta 

C 

176    IF(HOME)  THEN 

WRITE  (• ,  •) 'Problom  in  honing  Scanivalvos' 

ELSE 

WRITE  (• ,•) 'Problwt  in  stopping  Scanivalvos' 

ENOIF 

IF  HER.EQ2)  WRITE  (•.O'Dota  tronsfor  failtd  in  J-Boord' 

IF  (IER.EO.3)  WRITE  (•.•)' Fai I uro  in  dotoching  J-Board' 

IF  (HOME)  GOTO  16 

GOTO  26 

END 
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SUBROUTINE  SETUP(CLK_ID.OEV_IO.RATE.MCHANL.NPT.NPTS) 

C 

C   Subroutine  written  by  Moj  Richard  R  Ry I •» .  NPS  Student.  13  May  67 

C  Updated    15  Jul  87 

C 

C  This  subroutine  is  intended  to  attach  the  internal  clock  and  A-0  board 

C   to  the  main  program  and  allow  them  to  be  used  for  data  acquisition. 

C 

C   This  subroutine  should  be  used  with  subroutines  GETOATA,  DETACH 


C- 
C 

c- 


INCLUOE      SYS$LI8RARY:L10SET.FOR' 


Declare  local  variables 


C- 
C 

c- 


CHARACTER»4  CLK.DEV 
CHARACTER* 1  ANS 
INTEGER  SYS.STAT 
INTEGER  DEV_ID,CLK_ID 
INTEGER  NCHAN 
INTEGER  NPT.NPTS 
REAL    RATE. TIME 


'Status  ret  by  LIO  cal Is 
' LIO  device  ID 
'Channel  | 


Input  desired  clock  ond  sample  parameters 


WRITE(«.1) 
READ(».2)NCHANL 
10      WRITE(».3) 

READ(».»)RATE 
IF(RATE.GT.600.00)THEN 

WRITE(».4) 

GOTO  16 
ELSE 

IF(RATE.LT.0.ei)WRITE(«,5) 

IF(RATE.LT.0.ei)GOTO  10 
END  IF 
20      WRITE(».6) 

READ(«,7)  NPT 
IF(NPT.GT.5000)THEN 

WRITE(..8) 

GOTO  20 
ELSE 

IF(NPT.LE.0)WRITE(«.9) 

IF(NPT  LE.0)GOTO  20 
END  IF 

NPT  -  NPT  ♦  1 
TIME  -  F L0AT( NPT- 1) /RATE 
WRITE(».11)TIME. TIME/60. 
NPTS  -  2  •  NPT 
CLK_DEV  -  "KZA©- 
WRITE(«.12) 
READ(«.13)ANS 
IF  (ANS.EO. 'y'  OR  ANS.EQ. 'Y')  THEN 

WRITE(«.14) 

READ(«.15)CLK_DEV 
END  IF 


C- 

C 

C- 


C- 

c 
c- 


Attach  the  KWV11-C  Clock  ond  set  up  for  mopped  I/O 


SYS.STAT  -  LIO$ATTACH  (CLK_ID.  CLK_DEV.  LI0$K_MAP) 

IF(  NOT   SYS.STAT)  CALL  L ll$SIGNAL (XVAL(SYS.STAT) ) 

SYS.STAT  -  LIO$SET_R  (CLK_ID.  LlO$K_CLK_RATE .  1 .RATE) 

IF(  NOT   SYS.STAT)  CALL  L I8$SIGNAL(XVAL(SYS  STAT)) 

SYS.STAT  -  LI0$SET_I  (CLK_1D,  LIO$K_TRIG.  1.  LlO$K_IMMEDIATE) 
IF(  NOT   SYS.STAT)  CALL  L I B$S 1GNAL (XVAL( SYS.STAT) ) 


Attach  the  ADV11-D  ond  set  up  for  mopped  I/O 


SYS.STAT  -  L10SATTACH  (DEV.ID.  ' AZA0 ' .  LlO$K_MAP) 
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IF (NOT.  SYS.STAT)  CALL  LIB*SIGNAL(%VAL(SYS_STAT)) 

C 

C  S«t  up  th«  AOV  for  synchronous  transfer 

C       Set  channel  (nchanl) 

C       Coin  of  1 

C       Trigger  on  LIO$READ  and  fill  buffer  as  fast  as  possible 

C 

SYS  STAT  -  LIOSSET  I  (DEV.ID.  LlO$K_AD_CHAN.  1, NCHANL.) 

IF(.NOT.  SYS.STAT)  CALL  LIB$SIGNAL(XVAL(SYS_STAT)) 

SYS.STAT  -  LIO$SET  I  (DEV.ID.  LIO$K_AD_GAIN,  1.  1.) 

IF(.NOT.  SYS_STAT)  CALL  LIB$SIGNAL(XVAL(SYS  STAT)) 

SYS.STAT  -  LIO$SET  I  (DEV  ID.  LI0$K_TRIG.2 . LIO$K_CLK  POINT .CLK.ID) 
IF(.NOT.  SYS.STAT)  CALL  LIB$SIGNAL(XVAL(SYS_STAT) ) 

C 

C       Format  Statements 

C 

1  FORMAT( 1  x ,' Enter  channel  number  to  somple:  (0,1,8.9)') 

2  FORMAT(I2) 

3  FORMAT(/1x, • Input  desired  Clock  Rote  (606  Hz  -  0.01  Hz):1/) 

4  FORMAT(1 x . ' Fol lowi ng  directions!  Input  too  high . '//20x , 'Try  again'/) 

5  FORMAT( 1 x , ' Too  low!  Lots  of  time  on  your  hands? '//20x ,' Try  again'/) 

6  FORMAT (/1 x , 'No.  of  data  points:  (any  number  less  than  5.001)'/) 

7  F0RMAT(I4) 

8  FORMAT(/1 x . 'Hey !  Follow  directions!  Too  high.  Try  again!'/) 

9  FORMAT (/1 x . 'Hey !  Last  input  is  nonsense!  Try  again!'/) 

11  FORMAT (/1 x ,' The  approximate  run  time  for  each  port  will  be  * . f 6 . 1 » 
k  '    seconds  (or  '  ,  f  5 . 2 . '  minutes)'/) 

12  FORMAT (/1x, 'Do  you  want  to  change  the  Clock  Device  Nome  from  [KZA0]? 
*         "  (y/n)  [DEF  -  no]'/) 

13  FORMAT(A1) 

14  FORMAT(1x, 'Enter  Clock  Device  Name:  (K2A0.KZB0) ' ) 

15  FORMAT (A4) 

C 

RETURN 
END 
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SUBROUTINE  GETDATA(DEV_ID.NPT.NPTS. VOLTAGES) 


C 
C 
C 
C 
C 
C 
C 

c 
c 


Subroutine  written  by  Moj  Richard  R  Ry 1 •• .  MPS  Student,  28  Apr  87 

Update  15  Jul  87 

This  subroutine  ie  intended  to  acquire  the  data  from  the  A-0  board 
with  the  specified  parameters  from  subroutine  SETUP.  It  also  will 
translate  the  A— 0  data  into  voltages  to  be  converted  to  physical 
quantities  in  the  main  program. 

This  subroutine  should  be  used  with  subroutines  SETUP. DETACH, OUTPUTD 


C- 
C 


INCLUDE  'SYS$LIBRARY:LIOSrr.FOR' 


Declare  local  variables 


INTEGER  SYS_STAT 
INTEGER  DEV.ID 
INTEGER  DATA_LENGTH 
INTEGER  NPT.NPTS 


'Status  ret  by  LIO  calls 

1  LIO  device  ID 

'Number  of  bytes  of  data  read 

'Number  of  points  to  sample 


C 

c- 


Declare  data  buffers 


INTEGERS  RAW_DATA(5000) 
REAL     VOLTAGES (5090) 


15000  point  row  dato  buffer 
15000  voltages  output 


C  Get  a  raw_data  buffer  of  'npt'  values 
C       Uses  LIO$READ  to  read  the  A/D 
C       Note  length  of  buffer  is  in  bytes 
C       As  is  the  returned  data  length 


C 
C- 


SYS   STAT  -   LIO$READ(DEV    ID.    RAW  DATA.    NPTS.    DATA. LENGTH.    ) 
IF(.NOT.    SYS.STAT)    CALL    LIB$SIGNAL(XVAL(SYS_STAT) ) 

Convert    the    raw   data    to   voltages    using    LSP$FORMAT_TRANSLATE_ADC 


CALL    LSP$FORMAT_TRANSLATE_ADC(RAW_DATA.    VOLTAGES.    NPT 


.    ) 


RETURN 
END 
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SUBROUTINE  DETACH (CLK_IO.OEV.ID) 

C — 

C    Subroutine  written  by  Moj  Richard  R  Rylee.  MPS  Student  26  Apr  87 

C  Update  16  jun  87 

C 

C    This  subroutine  is  intended  to  detach  the  internal  clock  and  A-0  board 

C    from  the  main  program  and  allow  them  to  be  used  for  others  application. 

C 

C       This  subroutine  should  be  used  with  subroutine  SETUP 

C 

INCLUDE  ■SYSSLIBRARYiLIOSET.FOR' 

c.«..».« ••• • 

C        Declare  local  variables 

c..» • ••••• 

INTEGER  SYS  STAT        iStotus  ret  by  LIO  calls 

INTEGER  DEV_ID.CLK_ID    ! LIO  device  ID 

C. • • 

C       Detach  from  the  Clock 

C. • 

SYS  STAT  -  LIOSDETACH  (CLK.ID.) 

IF(.NOT.  SYS.STAT)  CALL  LIB$SIGNAL(XVAL(SYS_STAT) ) 

c...«« • 

C       Detach  from  the  A/D 

C.»« 

SYS.STAT   =    LIO$DETACH(DEV_ID.    ) 

IF(.NOT.    SYS.STAT)    CALL    LIB$SIGNAL(%VAL(SYS_STAT) ) 

WRITE(»,1) 
1  FORMAT (/1x. "Both   clock   and   A-D   device   are   detached'/) 

RETURN 

END 
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c- 
c 
c 
c 

c 

c 

c 

c- 

c 

c- 


c- 
c 

c- 


SU6R0UTINE  OUTPUTO(RATE.NPT.MCHANL.AilEAN.-STO_DeV.  J  .VOLTAGES) 

Subroutine  written  by  Uoj  Richard  R  Ry I •• .  NPS  Student  28  Apr  87 

Update  24  Jul  87 

This  subroutine  is  intended  to  output  diagnostic  results  to  a  file. 

This  subroutine  should  be  used  with  subroutine  SETUP.  DETACH.  GETOATA 

Declore  locol  variables 


INTEGER  NCHANL 

INTEGER  J 

INTEGER  NPT 

REAL»8  AMEAN.STD_OEV 

REAL  RATE. VOLTAGES (50Oe) 


'Channel  / 

{General  loop  indexes 

'Number  of  points 


Output  data  to  a  file 


OPEN(UNIT  -  16.  FILE  -  'SCANIDIAG  DAT '  .  STATUS  -  'N9') 

WRITE(ie.l)  J.  RATE.  NPT-1 .  NCHANI.  AMEAN.  STD.DEV 
1       FORMAT(/1x, 'Scan i valve  #'.  \l/"\  x,  'Clock  Rate  -' . f 6. 1 ,5x . 'Data  Points'. 
I         'Sampled  ■', i 4/1 x . 'Chonne I ' ,5x . "Mean ', 4x , 'Standard  Deviation'/ 
k  3x. i3.2f 12.4/) 

WRITE(10.») 'Scanival ve  f'.J 

DO  20  1-1 .NPT 
26  WRITE(ie..)'I*' .I.VOLTAGES(I) 

CLOSE(10) 

RETURN 

END 
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c« 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 

c 


SUBROUTINE  HWWA(NSAMPLES, ANSWER. BLUE_GREEN) 


PURPOSE: 

Dr  i  ves 


the  acquisition  of  the  LDV  data. 


ERROR  HANDLING: 

I ER  -  1  means  no  "fatal"  error. 

EXTERNAL  REFERENCES: 

NAME       DESCRIPTION  AND  SOURCE 

DVWA      Calls  DRV11-WA  device  driver 


to  acqui  re  A/D  data 


DEVELOPMENT 
DATE 
07/26/83 
12/12/85 

02/25/87 
05/06/87 
06/08/87 
07/15/87 

07/24/87 
08/29/87 


HISTORY: 
INITIALS 
TML 
CLH 

GBG 
CLH 
RRR 
RRR 


RRR 
RRR 


DESCRIPTION 
Adapted  from  DAPNT 

Modified  to  print  results  in  either 
octal  or  integer. 

Modified  for  use  with  DRV11-WA  and  VMS 
Modified  to  accept  "new"  DAMUX1 . 
Modified  to  subroutine  +  deleted  CLH  5/87  Mod 
Modified  to  calculate  velocities  +  optimized 
for  speed  of  execution 
Modified  to  view  Histogram 
Modified  to  occept  calibration  constant 


AUTHOR(S) :  Ted  Li chtenste i n ,  Informatics  General  Corp. 


IMPLICIT  NONE 

PARAMETER 
PARAMETER 
PARAMETER 

INTEGERS 

INTEGERS 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

LOGICAL 

CHARACTER*  19 

CHARACTER* 1 


NRAWMX-80000 

LUNTI-5 

LUNOUT-6 

IRAWDT(NRAWMX) 

1 1RAWDT( NRAWMX) . I 2RAWDT( NRAWMX) 

I ER . H I STO . Nf r . ONE , TWO 

IYESNO,  NRAWIN.  I.  BLUE.GREEN 

NCHANS.  NSAMPLES,  INDEX.  JNDEX 

NS . ANSWER . COUNT . END 

TIME (NRAWMX) . FR EQ ( NRAWMX ) .VELOC(NRAWMX) , I 3RAWDT(NRAWMX) 

MIX(2) .GAM(2) ,DEL(2) ,KAP(2) .0(2) , I (2) ,MEAN(2) . FRE0UENCY(2 .8001  ) 

FREQONE.FREQTWO,FRE(NRAWMX),FREQ1_AVE.FREQ2_AVE.CAL_CONST(2) 

DATSELCT    (3) 

DIRECTION 

ANS 


DATA  MIX(1).MIX(2)   /  -1.E-6.-1.E-6/ 
INCLUDE    '($SYSSRVNAM) ' 


100      CONTINUE 


C 

c 
c 
c 
c 


Fill  the  raw  data  buffer  with  this  bit  pattern: 

"1010101010101010" 

This  equals  125252(8).  -21846(10). 


200 


DO  200  1=1.  NRAWMX 

IRAWDT(I)  =  '  125252'0 
HRAWDT(I)  =  '  006000' o 
I2RAWDT(I)  =  ' 000000  o 
I 3RAWDT ( I )  =  0.0 

CONTINUE 


C 

c 

c 


Input    number    of    satiples    to   be   acquired 
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NCHANS  -  4 
3ee   WRITE  (LUNTI .  1B) 

Id     FORMAT  ('  Enter  number  of  samples:') 
READ  (LUNTI.  •)  NSAMPLES 
NRAWIN  -  NSAMPLES  •  NCHANS 

IF  (NRAWIN  .GT.  32767  .OR.  NRAWIN  .LE.  0)  THEN 
WRITE  (LUNTI ,  20) 
20       FORMAT  ('  The  number  of  samples  •  number  of  channels  must  be  32,767 
k  '  or  I  ess . ' ) 

GO  TO  300 
END  IF 

C  •  •  • 

C     Set  up  the  LDV  parameters  prior  to  acquisition 

C  ...... 

IF(BLUE_GREEN.EQ. 1)THEN 
GAM(1)  =  0.488 
GAM(2)  =  0.5145 
IF(CAL_CONST(1).E0.1.)  THEN 

WRITE(» ,*)' Enter  Blue  calibration  constant' 
READ(..«)  CAL_CONST(1) 
END  IF 
IF(CAL_C0NST(2).E0.1 .)  THEN 

WRITE(* .») "Enter  Green  calibration  constant' 
READ(..«)  CAL_C0NST(2) 
END  IF 
ELSE 

GAM(1)  «  0.5145 
GAM(2)  =  0.488 
IF(CAL_CONST(1).EQ.1 .)  THEN 

WRITE(* ,*) 'Enter  Green  calibration  constant' 
READ(..»)  CAL_CONST(1) 
END  IF 
IF(CAL_CONST(2).EQ.1 .)  THEN 

WRITE( ..»)' Enter  Blue  calibration  constant' 
READ(..«)  CAL_CONST(2) 
END  IF 
END  IF 
IF(MIX(1)  .LT  0.0)  THEN 

IF(BLUE_GREEN.EQ. 1 )WRITE(*  .»)' Input  BLUE  mixing  frequency  (MHz)' 
IF(BLUE  GREEN. EQ. 0)WRITE( •.•)' Input  GREEN  mixing  frequency  (MHz)' 
READ(..«)  MIX(1) 
END  IF 
IF(MIX(2)  LT  0  0)  THEN 

IF(BLUE_GREEN.EQ. 1 )WRITE(» .»)" Input  GREEN  mixing  frequency  (MHz)' 
IF(BLUE_GREEN.EQ.0)WRITE(»..),Input  BLUE  mixing  frequency  (MHz)' 
READ(..«)  MIX(2) 
END  IF 
I F( ANSWER  EO  0)  THEN 

WRI TE( • . • ) 'Do  you  wont  probe  volume  specifications?  (y/n)' 
READ(«. '(a) ' )  ANS 

IF(ANS  EO   Y'  OR.ANS.EQ. 'y' )  ANSWER  -  1 
END  IF 

WRITE(.,«)'  ' 
C 

C       Select  the  dota  to  be  returned. 
C 

DATSELCT  (1)  -  FALSE. 
DATSELCT  (2)  -  FALSE. 
DATSELCT  (3)  -   TRUE. 

C 

C      LDV  Dato  Acqui s 1 1  i on 

C 

ONE  -  1 

TWO  -  2 

OPEN (UN  I T-5. FILE- 'VELOCITY. DAT'  . STATUS- ' NEW ) 

DO  INDEX  -  1 .  NSAMPLES 

CALL  DAMUX  (NCHANS.  ONE.  DATSELCT,  ,  ,  I RAWDT .  I ER ) 
IF(IER   NE.  1)  WRITE(..*)'  IER  in  DAMUX  -'.IER 
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c  • 

C    Conversion  to  Frequency 
C  •  •  • •••< 


DO  1-1 .2 

IRAWDT(I)  -  INOT(IRAWDT(I)) 

I3RAWDT(I)  -  0.0 

I2RAWDT(I)  -  0 

HRAWDT(I)  -  0 

HRAWDT(I)  -  IIBITS(IRAWDT(I).10.4) 

I2RAWDT(I)  =  IIBITS(IRAK»DT(I),e.ie) 

I3RAWDT(I)  -  (2.»»I1RAWDT(I))  .(  f I  oat ( I2RAWDT( I ) ) ) 

TIME(I)  -  I3RAWDT(I)/32. 

FREQ(I)  -  ieee./TIME(I) 
END  DO 


FREQUENCY(1 .INDEX)  -  FREQ(1) 


FREQUENCY (2, INDEX)  -  FREQ(2) 
C  •  • 

C     Convert  frequency  data  to  velocities 


CALL  VELOCITY  (FREO.VELOC , TWO. MIX .ANSWER .GAM.DEL.KAP .0 . I ,Nf r .CAL.CONST) 
WRITE(5.30)  VELOC(1),VELOC(2) 
30      FORMAT(1X.2F11 .4) 
END  DO 
CLOSE(5) 

c...  ................................... 

C     Calculate  average  counter  frequency  and  print  to  sceen  for  diagnostics 
C  •  »  •  •  •  • • 

FREOONE  -00 
FREQTWO  -  0.0 
DO  34  I  -  1 .NSAMPLES 

FREOONE  -  FREOONE  +  FREQUENCY( 1 . I ) 
34       FREQTWO  *  FREQTWO  +  FREQUENCY(2 . I ) 
FREQ1_AVE  -  FREQONE/FLOAT(NSAMPLES) 
FREQ2.AVE  -  FREQTWO/FLOAT(NSAMPLES) 

WRITE(«.»)'  ' 

WRITE(«.») 'Average  counter  frequency  #  1  «',FREQ1_AVE 
WRITE(«,») 'Average  counter  frequency  #  2  = ' . FREQ2.AVE 
C  •  •  •  •  •  • • 

C     Write  Probe  Volume  dimensions  and  parameters  to  a  file 
C  • 


IF(ANSWER.EQ. 1)THEN 

OPEN(UNIT=15.FILE='Probe.dat ' . STATUS* ' NEW ' ) 

WRITE(15,«) 'Probe  Volume  Data' 

WRITE(15,»)' 


DO  40  I  -  1.2 

IF(BLUE_GREEN.EQ.1)  DIRECTION-  'Horizontal  Velocity' 
IF(BLUE  GREEN. EQ.0)  DIRECTION-  'Vertical  Velocity 
IF(I.EQ.2)  THEN 

IF(BLUE_GREEN.EQ.1)  DIRECTION-  'Vertical  Velocity   ' 
IF(BLUE_GREEN.EQ.0)  DIRECTION  -  'Horizontal  Velocity' 
END  IF 

WRITE(15,1)GAM(I) .DIRECTION , DEL( I ) ,KAP( I ) .D( I ) . I ( I )• . 001 ,Nf r 
FORMAT(/1X, 'Wave  length  of  light  <microns>  -  '.f7.4.'    (Measuring' 
4  a19 . ' ) '//1 x , ' Fr i nge  Spacing  <microns>  =  ' , f  10 . 2//1 x , ' Lens ' , 

k  '  angle  <deg>  «' . f 7 . 2//1 x , 'Probe  volume  diameter  <microns> 

k  f 10. 5//1 x ,' Length  of  measuring  volume  <mm>  -  ',f10.4//1x, 

4  Number  of  Fringes  =',i3/) 


WRITE(15.»)'- 


40        CONTINUE 
CL0SE(15) 
END  IF 

C  •  •  •  • •  i 

C     Option  to  print  a  histogram  of  velocity  data  on  screen  for  diagnostics 
<;••••••••• 

WRITE(«.«)'  ' 

WRITE(«.») 'Do  you  want  a  Histogram?  (1/0)' 
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READ(».»)   HISTO 

END  -   NSAJyFLES 

IF(HISTO.EQ.l)CALl  HISTOGRAM (COUNT . END. BLUE  GREEN) 


C     Again? 
C • 


WRITE(LUNTI,.)'  ' 

WRITE  (LUNTI.«)  'Acquire  LDV  Data  again?   (©-NO.  1-YES)  ' 

READ   (LUNTI,.)  IYESNO 

WRITE(».*)*  ' 

IF  (  IYESNO. EQ. 1  )  GO  TO  100 


RETURN 
END 
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c- 

c 

c 

c 

c 

c 

c 

c- 


SUBROUTINE  VELOCITY(FREO.VELOC.N.yiX,ANSWD».GAII.DEL.KAP.D«.  I  .Nfr. 
k  CAL.CONST) 

Subroutine  written  by  Moj  Richard  R  Rylee.  NPS  Student.  9  Jun  87 

Lost  Update  29  Aug  87 

This  subroutine  was  written  to  convert  the  LDV  raw  data  from  frequency 
to  ve I oc  i  t  ies. 

REAL    THETA. GAMMA. DELTA. d. f . E. KAPPA, Dsubm. Im.De . FREO(N) ,VELOC(N) ,GAM(2) 

REAL    SHIFT. MlX(2).DEL(2).CAL_CONST(2).KAP(2).Dm(2). 1(2) 

INTEGER  I ANS. ANSWER. Nfr 

IF(CAL_C0NST(1).EQ  0.0)  CAL.CONST(l)  «  1.6 

IF(CAL_CONST(2).EQ.0.0)  CAL_C0NST(2)  -  1.0 

PI  -  ASIN(1.0)«2. 

d  ■  13.  '  Bean  spac  i  ng 

f  -  600.  !  Focal  length 

E  «  1 .  !  Beam 

KAPPA  -  ATAND(E»d/2./f)        !  Lens 

DO  5  I  -  1  .2 

GAMMA  -  GAM(I) 

DELTA  -  CAL_CONST(I) 

De  -  1.1 

Dsubm  »  (4.»GAMMA»f/(PI»De»E))  ! 

Im  -  Dsubm  /  TAND(KAPPA)        ! 


(mm) 

of  transmitting  lens  (mm) 


expons i on  rat  i  o 
angle  (deg) 


Nfr  »  1 
DEL(I) 
KAP(I) 
Dm(I)  - 

l(I)  - 

SHIFT  = 

FREQ(I)  > 
VELCC(i) 

RETURN 

ENO 


/  De 


.27  •  d 

■  DELTA 

*   KAPPA  ! 

Dsubm  ! 

Im  ! 

FREO(I)  -  SHIFT  +  MIX(I) 
.  DELTA  •  FREO( i ) 


!  Wavelength  of  Loser  (microns) 
•  GAMMA  /  (2.«SIND(KAPPA))   !  Fringe  spacing 
!  Diameter  of  the  incident  laser  beam   (im 
Probe  volume  diameter  (microns) 
Probe  volume  length  (microns) 
Number  of  fringes 
Variobles  passed  to  output  file 


-  Bragg  shift  (MHz) 
MIX  -  Mixing  Freq  (Mr 
FREQ  -  Adjusted  Dopp 
(MHz) 


2) 

er 


Freq 
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c- 

c 

c 

c 

c 

c 

c 

c- 


SUBROUTINE  H I  STOGRAM(  COUNTER.  END.  BLUC_G8EEN) 


Subroutine  written  by  Moj  Richord  R  Ry i •• .  >e?S  Student.  16  Jun  87 

Updote  29  Aug  87 

This  subroutine  it  intended  to  creote  two  histogram*  of  velocity 
data  fro*  a  two  component  LOV. 


CHARACTER**     STRING 

CHARACTER*22  LABELX  t    variable*   for    the  oxi* 

CHARACTER* 17  LABELY  I    lab  lee 

CHARACTER. 29  LABELM 

INTEGER  STATUS. 1. 1 SHADE.  COUNTER.  D€> 

INTEGER  N(3).BLUE_GREEN 

REAL  XCONTROL(4).XMIN.XMAX 

REAL  YCONTROL(4).COUNT(8) 

REAL  XLOW(8) 

REAL  XHIGH(8) 

REAL  Y(8).VELOCITY(8288.2).X(8286) 

C 

N(1)   -   1  I    number   of   data   point* 

N(2)   -    1  I    number  of   columns  of   dato 

N(3)   -    1  I    number   of    rows  of    data 

C   ........ ••• 

C     Read  Velocity  Data  File 

C ••••• • 

OPEN ( UN I T-5. FILE-' VELOCITY. DAT '.STATUS- "OLD') 

IF(8LUE_GREEN. EO. 1 )THEN 
DO    10    I    -    1 .ENO 
16  READ(5.1)   VEL0CITY(I.1).VELCCITY(I.2) 

ELSE 

DO  28   1-1. END 
28  READ(5.1)   VELOCITY(I .2) .VELOCITY( I . 1 ) 

EX)IF 

C seeseeee.es 

C            Find    the   Min   and  Max   values 
C •• 

CLOSE(5) 

DO  88  J   -   1  . 2 

DO  38  K  -  1 , end 
38  X(K)  -  VELOCITY  (K.J) 

XMIN  -  1886.8 

XMAX  -  -1866.6 

DO  48  I  -  Lend 

IF(XMIN.GT.xm)  XMIN  -  XM) 
46  IF(XMAX.LT.X(I))  XMAX  -  X(I) 

C  ..........  • 

C     Initialize  X-Ax i s  controls 

c  .................  . 

XCONTROlM)   -  5.5  !     length   of    x-axis.    in    inches 

XCONTROL(2)   -  XMIN  !     left   most    value   on    x-axis 

XCONTROL(3)    -   XMAX  !     right    most    value    on    x-axis 

XCONTROL(4)    -    (XMAX-XMIN)/8 .  I     increment    of    plotted    points 

Ceeeeeesseee 

C  Find   y   data   from   sample   array 

C   ......... 

DO  56    I    -   1.8 

Y(I)   -  6.6 

XLOWf(I)  -  XMIN  +  FL0AT(I-1).XC0NTROL(4) 
IF(I.E08)THEN 
XHIGH(8)  -  XMAX 
ELSE 

XHIGH(I)-  XMIN  +  FL0AT(I)«XC0NTROL(4) 
ENDIF 
56  COUNT(I)  -  XMIN  ♦  FLOAT( I ).XC0NTROL(4) 

C 

XTOTAL  -88 
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c 
c 
c 


80 
90 


C 
C 
C 
1 
2 


OOM  1  -1,md 

XTOTAL  -  X(I)  ♦  XTOTAL 

IF(xm.GE.XMIN.AND.X(I).LT.COUNT(1))  Y(1)  -  Y(1)  ♦  1  e 
IF(X(l).CE.COUNTfn.AND.xn).LT.COUNT(2))  Y(2)  -  Y(2)  ♦  1.6 
IF(X(l).GE.C<XJNT(2).AND.X(l).LT.COUNTf3))  Yfj)  -  Y(3)  ♦  1.8 
IF(Xm.GE.COUNT(3)  AND.x(l ) . LT.COUNT(4))  Yf4)  -  Y(4)  ♦  1.6 
IF(X(l).GE.COUNT(4).AN0.Xfl).LT.C0UNT(5))  Yf5)  -  Y(5)  ♦  1.6 
IF(xn).GE.COUNT(5).AND.Xm.LT.COUNT(6))  Y(6)  -  Y(6)  +  1.6 
IF(xm.GE.COUNT(6).AND.X(l).LT.C<XJNTm)  Y(7)  -  Y(7)  ♦  1.6 
IF(X(I).GE.COUNT(7).ANO.XU).LT.COUNT(8))  Y(8)  -  Y(8)  +  1.6 
66         IF(X(1)  GE.COUNT(8).AN0.X(I).LT.XMAX)  Y(8)  -  Y(8)  ♦  1.6 

XMEAN  -  XTOTAL/FLOAT(END) 
C 

YMAX  -  6.6 
DO  76  I  -  1.8 

Y(D  -  Y(I)/FLOAT(END) 
IF(YMAX.LE.Y(I))YMAX  -  Y(I) 
76         Y(I)  -  166.  •  Y(I) 

WRITE(«.2)YVIAX.166.  , XMEAN 
c. ••••••••••••••••• 

C  Initial i  z«  Y  axis   controls 

C   .................    . 

YCONTROLM)  -  4.5  I  height  of  graph  in  inches 

YCONTROL(2)  -  6.8  I  lower  boundary  value 

YCONTROL(3)  -  YMAX* 126.  !  upper  boundary  value 

YCONTROL(4)  -  5.6  I  increment  of  plotted  points 

C  ........ •  •  •  • 

C     Histogram  labels 

C eeeeeeeeeeeee 

LABELX  -  '   VELOCITIES  (mps)    '  !  x-oxis  label 

IF (COUNTER. GT.6)LABELX  -  'FINAL  VELOCITIES  (mps)'  I  x-axis  label 

LABELY  -  'PERCENTAGE  OF  PTS*  !  y-oxis  label 

LABELM  -  'HORIZONTAL  VELOCITY  HISTOGRAM'        !  upper  label 
IF(J.EQ.2)  LABELM  -  '  VERTICAL  VELOCITY  HISTOGRAM- 
STRING  -  'IXSY'  !  mode  string 
I  SHADE  -  2                                   !  filling  of  graph 

C  .............*...*.*  . 

C     Set  up  coordinate  system  with  values 
C •••• 

IF(XMIN.GE.XMAX)GOTO  96 

CALL  LGPtPLOT (1. STRING.. 6. N. LABELX. LABELY. STATUS... 

XCONTROL . YCONTROL .. LABELM) 
IF  (STATUS  .NE.  LGP$_SUCCESS)  THEN 
CALL  LIB$SIGNAL(XVAL(STATUS)) 
ENDIF 


Call  routine  and  plot  histogram 

CALL  LGP$HIST(1.XLOW.XHIGH.Y. 8. STATUS.. ISHADE) 
IF  (STATUS  NE.  LGP$_SUCCESS)  THEN 
CALL  LIB$SIGNAL(XVAL(STATUS)) 
ENDIF 

PAUSE  'Waiting  for  viewing/printing  Histogram;  Type  CONT  when  finished' 
CALL  LGPSTERMINATE_PLOT(1)  I  terminate  plotting 

CONTINUE 
RETURN 

WRITER. •)'  ' 

WRITE(* ,•) '•••  No  Histogram  created: 
GOTO  86 


xmin>>xhigh  •••' 


Formats 


FORMATMX.2F11  .4) 
FORMAT(/1X.  'MAXII 


END 


F9.2) 


MUM  PERCENTAGE  OF  SECTORS  -  ' ,F6 . 1//1X. 'SAMPLE  MEAN  -  '. 
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SUBROUTINE  REDUCE_DATA( NSAMP, COUNTER. END, MEAN. STD_DEV1 .STD  DEV2. 
k  BLUE_GREEN) 

C 

C     Subroutine  written  by  Moj  Richard  R  Ryles.  NPS  Student.  12  Aug  87 

C  Update   29  Aug  87 

C 

C     This  subroutine  is  excludes  all  velocity  data  from  the  LDV  that  fall 

C     outside  two  standard  deviations  front  the  mean. 

C — 

REAL«8  STD.DEV1  ,STD_DEV2  .NEW_STD(2)  ,T0TAI_(2)  ,SUM(2)  .X.Y 

REAL  VELOC( 2.10060) .MEAN(2) .UP. DOWN. TOLERANCE 

REAL  NEW(2) . VEL0(2 . 20000) 

INTEGER  NSAMP.COUNT(2) . TOT , END.BLUE_GREEN 

C • 

C     Read  velocity  data  from  file 

C.  .......  ....................... 

OPEN ( UN I T=3, F I LE=' VELOCITY. DAT* . STATUS-' OLD ' ) 
DO   1 0    I    =    1 , NSAMP 
10  READ(3.1)VELOC(1 . I ) . VELOC(2 . I ) 

CLOSE(3) 

C 

C     Calculate  ronge  of  acceptable  values  i.e.  all  velocity 

C     data  that  falls  within  two  standard  deviations  of  the  mean 

C 

DO  30  J  -  1.2 

IF(J.EQ.1)THEN 

TOLERANCE  -  STD_DEV1 »2 . 
UP  -  MEAN(1)  ♦  TOLERANCE 
DOWN  -  MEAN(1)  -  TOLERANCE 
END  IF 
IF(J.EQ.2)THEN 

TOLERANCE  *  STD_DEV2«2. 
UP  -  MEAN(2)  +  TOLERANCE 
DOWN  -  MEAN(2)  -  TOLERANCE 
END  IF 

C 

C      Flag  each  point  outside  of  two  standard  deviations 

C  •  ........................  ...... 

DO  20  I  -  1 .NSAMP 
20  IF(VELOC(J. I).LT.DOWN.OR  VELOC(J. I ) .GT .UP)VELOC( J . I)«1000.0 

30    CONTINUE 

C  •  •  •  • •• 

C     Calculate  new  mean  and  standard  deviations 

C  •  •  • 

DO  60  J  -  1  . 2 

TOTAL(J)  -  0.D0 
COUNT(J)  «  0 
SUM(J)  -  0.D0 
DO  50  I  -  1  .NSAMP 

IF(VELOC(1 . I)  NE. 1000.0  AND  VELOC( 2. I ) . NE . 1000  0)  THEN 
TOTAL(J)  -  TOTAL(J)  +  VELOC(J.I) 
COUNT(J)  -  COUNT(J)  +  1 
END  IF 
50        CONTINUE 

1F(COUNT(J).EO0)  COUNT(J)  -  1 
NEW(J)  -  TOTAL(J)/FLOAT(COUNT(J)) 
DO  55  K  -  1 .NSAMP 

IF(VELOC(1 .K)  NE  1000.0  AND. VELOC( 2. K) . NE  1000  0)  THEN 
X  -  VELOC(J.K)  -  NEW(J) 
Y  -  X.X.(1 ./FLOAT (NSAMP-1 )) 
SUM(J)  -  SUM(J)  +  Y 
END  IF 
55       CONTINUE 

NEW_STD(J)  -  SORT(SUM(J)) 
60     CONTINUE 

C • • 

C     Output  new  volues  to  screen 

C • • 
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WRITE(«.»)'  ' 
tITE(». •)'*•< 


WRITE(» ,•) 'New  cleaned-up  data' 

I F ( BLUE_GREEN . EQ . 1 ) THEN 

WRITE(»,»)"U  di  r * ,NEW( 1 ) ,NEW_STD( 1 ) 
WRITE(«,»)'V  dir' ,NEW(2).NEW_STD(2) 

ELSE 

WRITE(».«)'U  di  r\NEW(2).NEW_STD(2) 
WRITE(v»)'V  di  r"  ,NEW( 1 ) ,NEW_STD( 1 ) 

END  IF 

MEAN(1 )  =  NEW(1) 

MEAN(2)  =  NEW(2) 

STD.DEV1  =  NEW_STD(1) 

STD_DEV2  =  NEW_STD(2) 

C  ...... 

C     Write  interim  data  file  containing  just  good  data 

c  ..........................  . 


OPENfUN^S.FILE^ INTERIM.DAT' , STATUS- ' NEW' ) 
DO  70  J  -  1  . 2 

DO  70  I  =  1  ,  NSAMP 

IF(VELOC(1 . I) .NE. 1000.0.AND.VELOC(2. I) .NE. 1000.0)  THEN 

WRITE(5.2)VELOC(J,I) 
END  IF 
70    CONTINUE 
CLOSE(5) 

C 

C     Write  final  velocity  dato  file  for  plotting 

(:••••••••••• 

OPEN( UN  I T=7 , F I LE= ' VELOC I TY . DAT '  . STATUS- ' NEW ' ) 
OPEN(UNIT«9. FILE-1 INTERIM.DAT1 , STATUS- 'OLD' ) 
TOT  -  COUNT(1)  +  COUNT(2) 
DO  80  I  =  1  ,  TOT 

IF(I .LE.COUNT(1))THEN 
READ(9.2)  VELO(I.I) 
ELSE 

READ(9.2)  VELO(2.I) 
END  IF 
80    CONTINUE 
CLOSE(9) 
END  =  COUNT(1) 
COUNTER  =  1 
DO  1 00  I  =  1 .  END 
100      WRITE(7,1)  VELO(1 , I ) . VELO(2 ,COUNT( 1 )+I ) 
CLOSE(7) 
RETURN 

C  ........ 

C     Format  statements 

c  .......  • 

1  FORMAT(1X.2F11 .4) 

2  FORMAT(1X.F11 .4) 
END 
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